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ABSTRACT 


k 


This  final  report  on  Contract  No.  F44620-73-C-0060, 
April  1,  1973  - March  15,  1978,  describes  the  work  car- 
ried out  under  the  terms  of  the  contract. 

A High  Gain  Long  Period  (HGLP)  seismograph  station 
was  installed  in  the  White  Pine  Copper  Mine  in  Ontonagan 
County,  Michigan.  The  station  went  through  a series  of 
stages  starting  with  a long  period  vertical  recording 
through  a phototube  amplifier  on  a strip  chart  and  ending 
as  an  essentially  standard  HGLP  station.  Because  of  the 
difficulties  of  servicing  and  recording  without  telemetry, 
the  station  was  closed  at  the  end  of  1976. 

Much  useful  information  regarding  the  engineering  and 
installation  of  HGLP  stations  in  mines  was  obtained.  The 
station  produced  many  months  of  useful  long  period  data. 

Concurrently  with  the  station  installation  two  ex- 
tensive research  efforts  in  long  period  seismology  were 
carried  on.  One  developed  a Rayleigh  wave  velocity  model 
for  ocean  basins  and  the  other  applied  long  period  P wave 
spectra  to  the  determination  of  seismic  moments. 


INTRODUCTION 


This  project,  lasting  nearly  five  years  from  April  1, 

1973  to  March  15,  1978,  involved  the  installation  of  a High 
Gain  Long  Period  (HGLP)  seismograph  station  in  the  White 
Pine  Copper  Mine  located  at  White  Pine,  Michigan  in 
Ontonagan  County.  The  station  was  developed  in  a series  of 
planned  steps.  It  was  initiated  in  1973  with  a vertical  com- 
ponent using  a phototube  amplifier  and  recording  on  a strip 
chart.  The  station  was  ultimately  expanded  to  three  com- 
ponents and  by  the  time  it  was  shut  down  in  1976,  it  was  in 
essence  a standard  HGLP  station. 

Various  innovations  were  made  in  the  installation. 

Perhaps  the  most  significant  was  the  manner  in  which  the 
horizontals  were  installed.  The  standard  large  steel  pres- 
sure containers  used  to  shield  the  horizontal  seismometers 
from  barometric  effects  were  modified  by  removing  their 
bottoms  and  cementing  the  rim  directly  to  the  sandstone  of 
the  mine  floor.  This  permitted  the  installation  of  the 
seismometer  directly  on  the  rock  rather  than  on  the  strained 
bottom  plate  of  the  container.  There  was  still  some  baro- 
metric shielding  but  as  the  horizontals  have  only  modest 
sensitivity  to  barometric  effects,  the  gain  in  stability, 
particularly  on  initial  set  up  or  at  times  of  adjustment, 
was  marked. 

The  history  of  the  station,  the  engineering  developments, 
and  the  various  details  of  the  installation  are  all  to  be 


found  in  the  semiannual  Technical  Reports.  Selected  sections 
of  them  are  given  again  in  this  final  report. 

Concurrentlyf  with  the  development  of  the  station  an  ex- 
tensive research  effort  was  carried  on  using  data  from  the 
HGLP  stations  and  from  long  period  seismograms  from  the  WWSSN 
stations . 


ASSOCIATED  RESEARCH 

The  research  in  long  period  seismology  led  to  two  Ph.D. 
theses.  One  by  Barbara  R.  Williams,  "Source  Parameters  for 
Large  Earthquakes  from  High  Gain  Long  Period  Body  Wave 
Spectra",  and  one  by  Frederick  J.  Mauk,  "A  Tectonic-Based 
Rayleigh  Wave  Group  Velocity  Model  for  Prediction  of  Dis- 
persion Character  through  Ocean  Basins". 

The  first  deals  with  the  use  of  long  period  spectra  of 
P and  PP  to  determine  the  moment  of  large  earthquakes.  A 
large  number  of  spectra  were  calculated  using  both  Fast 
Fourier  Transform  (FFT)  and  Maximum  Entropy  Spectra  Analysis 
(MESA)  methods.  It  became  evident  that  while  an  estimate  of 
the  moment  could  be  obtained,  the  noise  (probably  both  in- 
strumental and  earth)  above  100  seconds  was  too  large  to  make 
this  a completely  satisfactory  method.  The  corner  frequency 
could  not  be  determined  reliably  but  the  determinations  were 
consistent  with  moments  derived  from  field  observations  and 
aftershock  distribution.  An  extensive  elaboration  of  Brune's 
AR  method  was  made  and  the  method  further  adapted  to  WWSSN 
records  and  oceanic  paths.  Moments  calculated  by  the  AR 
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method  were  also  used  for  moment  comparison. 


The  second  thesis  developed  a new  model  to  predict 
Rayleigh  wave  group  velocity  dispersion  (periods  30  through 
100  seconds)  through  ocean  basins.  It  is  based  on  a 5®  x 5° 
integration  of  weighted  Rayleigh  group  velocities  for  the 
2592  geographic  grid  elements  constituting  the  earth's  sur- 
face. 

Both  theses  are  available  from  University  Microfilms, 

Ann  Arbor,  Michigan.  Appendix  I is  a draft  manuscript  from 
Dr.  Williams'  thesis  and  Appendix  II  is  a manuscript  sub- 
mitted to  JGR  from  Dr.  Mauk's  thesis.  Appendix  III  includes 
various  abstracts  presented  at  scientific  meetings. 

A detailed  log  of  all  events  recorded  at  the  White  Pine 
station  (WPM)  has  been  made  and  is  on  file  at  the  Seismo- 
, logical  Observatory,  Department  of  Geology  and  Mineralogy, 

('■.  University  of  Michigan.  Appendix  IV  is  a Catalog  of  Events 

for  the  period  July  1,  1974  through  April  30,  1975. 

DEVELOPMENT  OF  THE  HGLP  STATION  WPM 

The  location  of  the  station  was  chosen  because  of  the 
availability  of  a suitable  site  in  the  White  Pine  Copper 
Mine  in  Ontonagan  County,  Michigan  at  a geographic  location 
such  that  relatively  uncomplicated  paths  were  thought  to  exist 
(and  indeed  do)  across  the  arctic  from  Eurasia. 


The  location  was  89  meters  below  the  surface  at  an  ele- 
vation of  193.5  meters  in  the  Parting  Shale  member  of  the 


Nonesuch  formation  of  preCambrian  age.  The  geographic 
coordinates  of  the  vault  are  45®  45'  05. "5  N,  89®  33'  17. "8  W. 
The  location  is  in  drift  S13,  W4  of  section  31-A  of  the  mine. 

A stratigraphic  section  and  map  of  the  vault  are  shown  in 
Appendix  IV,  pages  5 and  6. 

This  portion  of  the  mine  is  a room  and  pillar  operation 
which  leaves  about  half  of  the  rock  to  support  the  overlying 
strata.  In  both  the  initial  area  for  the  vertical  (vault  1) 
and  the  expanded  area  (vault  2)  to  accommodate  the  horizontals, 
supporting  pillars  were  used  as  part  of  the  enclosure  and  the 
open  areas  closed  with  plywood  sealed  with  six  inches  of  poly- 
urethane foam. 

Within  the  enclosed  area  containing  the  seismometers,  the 
temperature  of  the  mine  fluctuated  only  a few  degrees  from 
55®F.  The  pressure  tanks  were  filled  with  spun  glass  insula- 
tion to  reduce  thermal  effects.  During  the  summer  months, 
the  humidity  was  high  (95%)  but  was  considerably  lower  in  the 
winter. 

The  following,  referring  to  the  initial  installation  of 

the  vertical,  is  excerpted  from  the  Technical  Report  for  the 

period  1 April  - 30  September  1973. 

Before  the  seismometer  pressure  tank  could  be  installed, 
the  floor  was  cleaned  and  all  loose  debris  removed.  Most 
of  the  floor  was  found  to  be  unsuitable  for  the  instal- 
lation of  the  one  meter  tank,  either  because  of  extensive 
cross  jointing  or  more  often  because  normal  production 
blasting  highly  fractures  the  floor. 
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To  reduce  the  problems  of  tilting  normally  associated 
with  the  curing  of  concrete,  the  tank  was  installed 
directly  on  bedrock  with  a minimal  amount  of  cement  to 
seal  the  tank  bottom  to  the  floor.  Four-threaded  cad- 
mium-plated rods  were  cemented  38  centimeters  into  bed- 
rock with  epoxy  resin  to  secure  the  tank  bottom.  Be- 
fore mounting  the  tank,  the  bottom  was  prestressed 
into  a dome  with  maximum  displacement  of  2 centimeters 
at  the  center.  The  tank  was  then  set  on  a slurry  of 
sand  and  cement  and  bolted  to  the  floor.  The  prestress- 
ing tension  was  then  gradually  released  over  the  next 
twenty-four  hours. 

In  addition  to  the  AC  electrical  system  provided  by  the 
mine,  data  and  control  cabling  had  to  be  installed.  To 
minimize  pressure  leakage  into  the  seismometer  tank,  all 
cables  entering  the  tank  were  stripped  of  insulation, 
sealed  with  General  Electric  Glyptal  1201  and  potted  in 
3-M  Scotchcast  resin  plugs.  These  resin  plugs  were  then 
seated  in  Marsh-Marine  fittings  to  insure  pressure  integ- 
rity of  the  sealed  system.  Control  cabling  consisted  of 
Belden  6-palr  shielded  instrument  cable.  All  other  wiring 
was  dene  with  individual  shielded-pair  cable  to  reduce 
cross-talk  between  pairs  and  maintain  integrity  of  the 
data  lines. 

The  seismometer  used  at  station  WPM  is  a Geotech  model  S-11 
which  has  been  modified  to  accommodate  a Sprengnether  model 
VCT-210  displacement  transducer.  The  natural  frequency  of 
the  instrument  has  been  set  at  0.003  Hertz  (T^=30  seconds). 

A preliminary  investigation  employing  a Geotech  model  5240A 
phototube  amplifier  was  run  from  23  August  to  25  September 
to  determine  the  feasibility  of  a permanent  installation. 

In  mid-October  the  phototube  amplification  system  will  be 
replaced  with  a Geotech  model  AS-650  low  noise  amplifier 
incorporating  telemetry  capabilities. 
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The  output  of  the  phototube  amplifier  was  recorded  at  two 
gain  settings  on  Esterline  Angus  model  S601-S  six-inch 
strip  chart  recorders.  Output  from  the  displacement  trans- 
ducer system  was  recorded  on  an  Esterline  Angus  model 
EllOl-S  ten-inch  strip  chart  recorder. 

Time  marks  are  provided  by  a Sprengnether  S-100  chronometer 
which  was  determined  to  have  a constant  drift  rate  of  -4000 
milliseconds  per  day.  All  power  for  amplifiers  and  re- 
corders is  controlled  by  a Sorensen  ACR-3000  power  regu- 
lator which  reduces  the  problems  of  power  fluctuations  at 
the  station.  The  most  severe  power  fluctuations  occur  when 
a mine  compressor  in  the  area  is  turned  on.  This  causes  a 
momentary  10  volt  power  dip  followed  by  a transient  surge. 

The  power  regulator  reduces  the  10  volt  dip  to  2 volts  and 
eliminates  the  associated  surge. 

Continuous  seismic  recordings  were  run  from  this  pre- 
liminary installation  through  the  period  August  23  to  Sep- 
tember 25,  1973.  Two  gain  levels  (40K  at  45  seconds  and 
400K  at  45  seconds)  were  used.  Figure  1 shows  the  response 
curve  at  400K  together  with  those  for  other  high  gain  stations. 
Figure  2 shows  the  seismogram  tor  a magnitude  5.3  event  from 
the  north  Atlantic  followed  by  a smaller  event  of  unknown 
location. 

More  details  and  illustrations  of  the  initial  installation 
can  be  found  in  the  first  Technical  Report. 

The  preliminary  installation  was  shut  down  from  September 
25,  1973  to  March  2,  1974,  while  awaiting  equipment  to  tele- 
meter the  signals  from  the  mine  to  the  surface.  Minimal  re- 
centering  was  required  to  return  to  operation.  A permanent 
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hard  wire  telemetry  link  was  installed  at  this  time  to  permit 
recording  at  one  of  the  mine  research  buildings  on  the  surface. 

In  the  summer  of  1974,  two  long  period  horizontals  were 
installed.  (See  Technical  Report  1 April  - 30  September,  1974)  . 
A second  vault  had  to  be  constructed.  The  following  is  from 
the  Technical  Report. 

All  loose  floor  debris  was  removed  before  wall  construc- 
tion was  begun.  This  substantially  reduced  the  difficulty 
of  selecting  tank  installation  sites  since  the  entire 
vault  2 floor  was  exposed.  Although  production  blasting 
had  reduced  the  usable  floor  space,  as  in  vault  1,  two 
unfractured,  unjointed  sections  of  the  Copper  Harbor  Con- 
glomerate were  deemed  suitable  for  tank  placement.  These 
sections  of  floor  were  meticulously  cleaned  and  any  loose 
fragments  removed.  This  resulted  in  relatively  smooth, 
impermeable  tank  foundations. 

Unlike  other  HGLP  installations,  where  prestressing  of 
the  tank  bottom  was  necessary  to  insure  proper  bonding  of 
the  metal  to  concrete  beneath  it,  we  chose  to  minimize 
tank  tilt  effects  on  the  horizontal  seismometers  by  de- 
coupling the  two  components.  Before  installation,  the 
tank  bottoms  were  cut  out  leaving  only  a rim  four  inches 
wide  for  bonding.  Four  cadmium-plated  threaded  rods  for 
each  tank  were  then  cemented  38  centimeters  into  bedrock 
using  an  epoxy  resin.  The  use  of  epoxy  bolts  rather  than 
mechanical  bolts  was  suggested  by  mine  personnel.  Both 
resin  and  mechanical  bolting  are  widely  employed  in  the 
White  Pine  Mine,  and  resin  bolting  has  been  shown  to  be 
more  effective  than  mechanical  bolting.  Since  it  was 
unnecessary  to  insure  a bottom  contact  with  concrete, 
prestressing  of  the  tank  was  unnecessary.  The  tank  was 
then  set  onto  a slurry  of  sand  and  cement  and  bolted 
securely  to  the  bedrock.  By  wobbling  the  tank  before 


final  bolting,  it  was  possible  to  reduce  the  average  thick- 
ness of  the  cement  to  less  than  5 millimeters.  The  cement 
was  then  permitted  to  cure  for  three  days  undisturbed. 
Following  the  curing  period,  the  inside  of  the  tanks,  bed- 
rock and  concrete  junctions  were  triple-sealed  with  epoxy. 

The  concrete  outside  the  tank  was  also  sealed.  Each  epoxy 
coat  was  permitted  to  harden  thoroughly  before  the  next 
was  applied.  Finally,  the  interior  and  exterior  of  each 
tank  was  coated  with  zinc  chromate  to  inhibit  rust  formation. 

All  but  two  ports  entering  the  tanks  were  sealed  with  teflon- 
coated  plugs  which  were  tightened  and  then  coated  internally 
and  externally  with  General  Electric  Glyptal  1201.  The  re- 
maining two  ports  were  fitted  with  Marsh-Marine  connectors 
which  clamped  onto  molded  wire  plugs. 

To  minimize  pressure  leakage  into  the  seismometer  tanks, 
all  cables  entering  the  tanks  were  stripped  of  insulation, 
sealed  with  General  Electric  Glyptal  1201,  and  potted  in 
3-M  Scotchcast  resin  plugs.  The  resin  plugs  were  then 
seated  in  the  Marsh-Marine  connectors  to  insure  pressure 
integrity  of  the  sealed  system.  Several  different  castings 
of  the  plugs  were  necessary  before  usable  plugs  could  be 
obtained.  It  was  found  that  the  3-M  Scotchcast  resin  became 
brittle  and  shattered  upon  removal  from  the  molds  if  per- 
mitted to  harden  in  the  ambient  mine  ci^nditions . Good  plugs 
could  be  obtained  only  when  the  molds  were  continuously 
warmed  with  heat  lamps  during  the  hardening  process.  Control 
cabling  consisted  of  Belden  6-palr  shielded  instrument  cable. 
All  other  wiring  was  done  with  individual-shielded  pair 
cable  to  reduce  cross-talk  between  pairs  and  maintain  in- 
tegrity of  the  data  lines. 

The  horizontal  seismometers  employed  at  station  WPM  are  Geo- 
tech  model  8700C's  which  have  been  modified  to  accommodate 
Sprengnether  model  VCT-210  displacement  transducers.  The 
natural  frequency  of  each  Instrument  has  been  set  at  0,033 
Hertz  (T^=30  sec).  Following  calibration,  each  seismometer 
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with  its  associated  displacement  transducer  electronics 
was  thoroughly  insulated  with  spun  glass  and  then  sealed 
in  its  tank.  With  no  physical  connection  between  the 
tank  and  seismometer,  C-clamps  could  be  tightened  until 
the  seismometers  "quieted  down".  There  was  a direct 
connection  between  the  output  noise  and  effectiveness  of 
the  tank  seal,  and  the  noise  level  would  dramatically  de- 
crease when  a good  seal  was  effected.  Since  there  was  no 
direct  connection  between  tank  and  seismometer,  stresses 
induced  in  the  tanks  during  clamping  were  not  transmitted 
as  tilts  to  the  seismometers.  Hence,  it  was  possible  to 
tighten  or  loosen  clamps  or  add  additional  clamps,  as 
necessary,  without  affecting  the  seismometers.  Through- 
out the  clamping  procedure  neither  seismometer  needed 
recenterlng. 

In  addition  to  the  seismometers,  a subsurface  control  area 
has  been  installed  in  the  northeast  corner  of  vault  1. 

Three  phototube  amplifiers  have  been  made  operational  as 
well  as  a Geotech  AS-650  solid  state  amplifier.  Recording 
of  horizontal  seismograms  was  done  on  two  Esterline  Angus 
model  S601-S  six-inch  strip  chart  recorders  in  the  sub- 
surface control  area  throughout  the  summer  of  1974.  In- 
strumentation in  the  subsurface  control  rack  consists  of 
a boom  position  monitor  and  reposition  motor  drive  panel, 
a calibration  panel,  a waveform  generator,  a Sprengnether 
S-100  chronometer,  and  power  regulation  and  distribution 
devices.  All  boom  monitoring  and  repositioning  for  both 
vertical  and  horizontal  seismometers,  as  well  as  cali- 
bration, must  be  effected  from  this  rack.  No  calibration 
or  repositioning  capability  exists,  or  is  likely  to  exist, 
at  the  surface  recording  site.  The  telemetry  link  to  the 
surface  recording  site  for  frequency  modulated  signals  is 
also  accessible  at  a panel  in  the  subsurface  control  area. 

As  of  October,  1974,  the  installation  of  a preliminary 
station  had  been  completed  and  it  was  possible  to  make  an 
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assessment  of  what  had  been  learned  and  to  feel  assured  that 


White  Pine  was  a more  than  suitable  site.  The  following 

lengthy  excerpt  from  the  Technical  Report  1 April  1974  - 

30  September  1974,  summarizes  this  information.  It  will  be 

noted  that  reference  is  made  to  a short  period  vertical  which 

was  operated  intermittently  as  a test  of  the  site. 

a)  Comparison  of  Geotech  AS-650  solid  state  amplifier 
and  phototube  amplifier  - In  June,  1974,  the  long  period 
vertical  instrument  was  switched  from  phototube  amplifier 
to  a Geotech  AS-650  solid  state  amplifier.  This  permitted 
recording  of  long  period  vertical  (in  addition  to  short 
period  vertical)  seismic  signals  on  the  surface  since  the 
output  of  the  AS-650  was  multiplexed  to  a 1020  Hertz 
carrier.  The  successful  recording  of  the  short  period 
vertical  signal  on  the  surface  has  been  discussed  pre- 
viously (see  Technical  Report  1 October  1973  - 31  March 
1974).  No  interference  or  direct  coupling  between  the 
long  period  and  short  period  systems  was  observable  on 
nearly  continuous  records  from  June  to  September,  1974. 

Further  discussion  of  the  short  period  system  is  given 
in  part  b)  of  this  section. 

The  Geotech  AS-650  amplifier,  unlike  the  phototube  ampli- 
fier, does  not  have  a 0,167  Hertz  notch  filter  incorpo- 
rated. This  resulted  in  high  six-second  microseismlc 
noise,  especially  during  "storm"  periods.  This  noise 
proved  to  be  the  principal  gain  limiting  factor  for  the 
WPM  long  period  vertical  system.  The  vertical  magnifi- 
cation ranged  from  88K  during  periods  of  intense  micro- 
seismic  activity  to  lOOK  during  quiet  "noise"  periods. 

The  bottom  trace  of  Figure  3 illustrates  the  signal 
appearance  during  an  intense  noise  period.  To  reduce  this 
effect,  a galvanometric  preconditioning  scheme  was  employ- 
ed (see  Pomeroy  et.  al.;  B.S.S.A.,  Vol.  50,  pp  135-151, 

1960).  A Leeds  and  Northrop  critically  damped  0.167  Hertz 
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Vertical  seismogram  Station  WM  To*=30  sec.  Geotech  AS-650  amplifier 
Gain  88,000  without  filter  galvanometer 
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galvanometer  was  inserted  in  series  between  the  seis- 
mometer and  the  AS-650  amplifier  (Pomeroy  P2F  configura- 
tion). The  galvanometer  in  this  mode  acted  as  a band 
rejection  filter,  preconditioning  the  signal  before 
amplification.  This  resulted  in  a six-second  micro- 
seismic  noise  reduction  by  a factor  of  two  (upper  trace. 

Figure  3).  The  use  of  the  galvanometric  band  rejection 
filter  scheme,  therefore,  is  being  employed  with  the 
AS-650  amplifier  at  WPM  on  HGLP  vertical  signals.  Two 
gain  levels,  lOK  and  lOOK,  are  being  recorded  for  the 
HGLP  vertical. 

During  the  period  1 through  4 August  1974,  a comparative 
test  was  run  using  the  vertical  seismometer.  The  signal 
from  one  data  coil  was  amplified  using  the  Geotech  AS-650 
solid  state  amplifier.  Simultaneously,  the  signal  from 
the  other  data  coil  was  amplified  using  the  Geotech  5240-A 
phototube  amplifier.  The  AS-650  employed  the  galvanometric 
band  rejection  filtering  technique  previously  described 
while  the  phototube  amplifier  employed  the  Geotech  6824-15 
special  filter.  Three  of  the  seismograms  recorded  in  this 
manner  are  included  in  this  report  (see  Figures  7-9).* 

This  comparison  revealed  several  important  facts:  most  of 
the  longer  period  (20  to  100  seconds)  noise  was  being 
generated  within  the  phototube  system  since  it  did  not 
appear  to  be  higher  for  periods  in  the  50  to  60-second 
range  from  the  phototube  amplifiers,  substantial  distortion 
of  10  to  20-second  periods  exists;  usually  body  phase 
arrivals  are  substantially  more  recognizable  on  the  AS-650 
records.  We  believe,  therefore,  that  the  seismograms  from 
the  AS-650  are  more  useful  to  seismologists  than  those 
from  the  Geotech  5240-A  phototube  amplifiers.  Hopefully, 
the  Ithaco  replacements  will  exhibit  more  of  the  character- 
istics of  the  AS-650  than  the  5240-A. 

*Flgure  7 is  shown  as  Figure  4 in  this  report. 
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b)  Suitability  of  the  White  Pine  Mine  for  a complete  HGLP 

i ■ 

Installation  - The  White  Pine  Copper  Company  experienced  a 
strike  shutdown  during  the  latter  two  weeks  of  August. 

During  this  period,  only  routine  maintenance  of  underground 
and  surface  facilities  was  undertaken  by  limited  mine 
personnel.  This  was,  therefore,  an  ideal  time  to  see  what 
affect  the  mining  operations  had  on  the  seismic  signals 
at  WPM.  There  was  no  noticeable  affect  on  any  of  the  three 
HGLP  instruments;  that  is,  there  appeared  to  be  no  reduction 
of  long  period  noise.  The  short  period  vertical  system 
underwent  substantial  signal  enhancement.  During  active 
mining  periods,  the  short  period  vertical  is  operable  at 
a maximum  gain  of  25K.  The  gain  limiting  factor  is  a per- 
sistent 4 Hertz  resonance  associated  with  mining  activities 
and  possibly  surface  milling  operations.  During  the  shut- 
down period,  the  gain  could  be  raised  by  at  least  a factor 
of  4 and  possibly  greater.  It  is,  therefore,  apparent 
that  to  run  a high  gain  vertical  short  period  instrument, 
it  must  be  remote  from  the  active  mining  area  of  the  White 
Pine  Copper  Company.  Prior  to  the  installation  of  the  HGLP 
vertical  instrument  in  1973,  a micro-earthquake  survey  was 
run.  Gains  in  excess  of  lOOK  could  be  obtained  at  a dis- 
tance of  10  miles  from  the  White  Pine  mining  operations. 

The  environmental  stability  of  the  WPK  ■seismic  vaults  has 
been  discussed  in  previous  technical  reports  (see  TR  1 April 
- 30  September  1973;  TR  1 October  1973  - 31  March  1974). 

At  that  time  the  long  period  vertical  instrument  had  been 
successfully  operated  at  gains  in  excess  of  400K.  The 
horizontal  instruments,  installed  during  June  and  July, 

1974,  also  have  proved  to  be  very  stable.  The  horizontal 
instruments,  in  fact,  have  been  far  more  stable  than  the 
vertical.  They  require  less  frequent  recentering  than  the 
vertical  and  have  gone  for  periods  of  a week  or  more, 
immediately  following  installation,  during  which  there  was 
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no  apparent  drift  at  all.  Presumably  this  stability  is  a 
reflection  of  the  decoupling  of  the  seismometer  from  its 
pressure  tank.  The  stability  of  the  horizontals  is  dis- 
rupted, however,  when  personnel  enter  vault  2 (possibly  a 
thermal  agitation  affect  or  actual  vault  tilt),  and  so 
entrance  is  restricted  to  major  maintenance  only.  Figure 
5b*  illustrates  the  response  curves  of  the  three  WPM  HGLP 
instruments.  Five  selected  events  recorded  at  WPM  are 
illustrated  in  Figures  10-14.** 

Difficulties  with  the  installation  during  this  report 
period  are  as  follows: 

1)  Stability  of  the  power  supplied  by  the  White  Pine 
Copper  Company  both  on  the  surface  and  underground  has 
become  suspect.  During  this  report  period,  the  under- 
ground installation  has  experienced  no  less  than  four  major 
power  outages,  the  surface  installation  more  than  ten. 

Damage  of  the  AS-650  amplifier  occurred  due  to  surging  when 
power  was  resumed.  To  prevent  recurrence  of  this,  an  ade- 
quate power  Isolation  system  is  required  both  above  and 
underground . 

2)  Recording  of  horizontal  seismograms  was  accomplished  in 
vault  1 of  the  underground  faciltiy.  To  complete  the  in- 
stallation, these  signals  must  be  telemetered  to  the  sur- 
face. This  includes  displacement  transducer  signals  as 
well  as  velocity  transducer  signals. 

3)  Many  problems  exist  with  the  Astrodata  data  logger 
system,  which  must  be  corrected  with  the  assistance  of  the 
Albuquerque  Seismological  Center  personnel. 

Following  the  completion  and  testing  of  the  "experimental" 
station  described  above,  it  was  decided  to  bring  the  station 

*Figure  5 in  this  report 

**Figures  10  and  13  are  shown  as  Figures  6 and  7 in  this  report. 
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Vertical  (high  and  low  gain),  North-South,  and  East-West  HCLP  records.  Event  at  Tadyblk-Sinklang 
Border  Region,  39.4  N,  73.9  E,  normal  depth  on  August  11,  1974  at  070208.5  UCT,  in>«5,8,  lla-5.4. 
For  HCLP  Kast-Weat  records,  advance  rate  Is  15  as  per  sinute  instead  of  13.5. 


into  a mode  as  nearly  like  the  standard  HGLP  stations  as 
feasible.  During  the  period  covered  by  the  Technical  Re- 
port 1 October  1974  - 31  March  1975,  planning  and  engineer- 
ing went  on  and  the  timing  system  was  upgraded.  Work  was 
initiated  with  the  Albuquerque  Seismological  Center  of  the 
U.  S.  Geological  Survey  through  which  a new  telemetry 
system,  Sprengnether  VCT-210  seismometer  systems,  and  Ithaco 
amplifiers  - all  standard  for  HGLP  stations  - were  readied 
for  installation.  An  uninterruptible  power  system  was  also 
readied  by  Albuquerque  to  avoid  the  unreliability  and  voltage 
irregularity  of  the  commercial  power  at  White  Pine. 

In  early  December  of  1974,  the  LPZ  velocity  signal  was 
transferred  from  strip  chart  to  rectilinear  helicorder  re- 
cording which  improved  the  continuity  of  recording. 

During  the  summer  of  1975,  much  of  the  work  to  convert 
the  station  was  completed.  The  following  excerpts  from  the 
Technical  Report  1 April  - 30  September  1975,  give  the  main 


details . 


Subsurface  instrumentation  - During  this  report  period 
the  following  alterations  of  the  White  Pine  subsurface 
installation  have  been  effected:  (1)  complete  rewiring 
of  the  vaults  with  shielded  solid  two-conductor  cable 
replacing  stranded  and  shielded  multiconductor  cables; 

(2)  installation  of  an  uninterruptible  power  system; 

(3)  replacement  of  the  phototube  and  AS-650  amplifiers 
with  Ithaco  solid  state  amplifiers;  (4)  installation  of 

a six-component  VCO  system  for  multiplexing  both  velocity 
and  displacement  data. 
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Rewiring  of  the  White  Pine  subsurface  vaults  was  accomp- 
lished in  July  at  the  suggestion  of  Albuquerque  Center 
personnel.  Difficulty  in  locating  the  source  of  high 
frequency  noise  induced  in  the  signals  prompted  this 
decision.  All  shielded  stranded  two-conductor  and 
shielded  multiconductor  cables  connecting  the  HGLP  tanks 
were  replaced  with  shielded  solid  two-conductor.  Although 
this  did  not  result  in  a significant  noise  reduction,  it 
considerably  improved  the  ease  of  subsurface  trouble  shoot- 
ing. 

A WWSSN-type  uninterruptible  power  system  consisting  of  a 
LaMarche  charger  and  Nicad  batteries  was  provided  by  the 
Albuquerque  Center  in  June.  Prior  to  this  time,  the  vault 
power  was  supplied  directly  by  available  circuits  in  the  mine. 

Replacement  of  the  Phototube  amplifiers  and  the  AS-650 
amplifier  with  Ithaco  KGLP  amplifiers  was  effected  in  July. 

With  this  alteration,  the  White  Pine  frequency  response  of  the 
velocity  outputs  will  match  all  other  HGLP  observatories. 

A six- component  VCO  system  with  active  mixers  was  added  to  the 
subsurface  instrumentation  in  June.  Initial  testing  of  the 
multiplexing-demultiplexing  portion  of  the  circuit  revealed 
excessive  high-frequency  noise  which  would  severely  limit  the 
gain  of  the  analog  records.  Although  inexperience  with  telemetry 
equipment  led  to  a random  logical  approach  to  the  noise  reduc- 
tion, eventually  all  major  sources  of  interference  were  located 
and  minimized.  Three  major  sources  of  Interference,  all 
associated  with  the  telemetry  circuits,,  were  responsible  for 
approximately  95%  of  the  noise.  Over-driving  the  telephone 
line  with  voltages  too  high  for  six  center  frequencies  was 
found  to  be  the  most  significant  noise  source.  When  the 
Individual  output  levels  were  adjusted  so  that  the  combined 
line  input  was  at  the  0 db  level,  the  noise  reduced  by  an 
order  of  magnitude.  The  second  noise  source  was  attributable 


-22- 


to  an  unfortunate  choice  of  center  frequencies.  Algebraic 
addition  of  several  carrier  frequency  pairs  was  found  to 
yield  resultants  which  Interfere  with  a third  channel. 

Avoiding  this  difficulty  using  constant  bandwidth  telemetry 
is  extremely  difficult  and  cannot  be  completely  eliminated 
but  only  minimized.  At  White  Pine  this  noise  reduction  was 
accomplished  by  splitting  the  six  components  into  two  mixed 
three-component  line  drivers  and  telemetering  the  signals 
on  two  separate  telemetry  lines.  In  addition,  the  center 
frequencies  were  adjusted  to  minimize  interchannel  inter- 
ference. The  third  noise  source  was  found  to  be  the  tele- 
metry line  itself . The  communication  line  provided  by  the 
White  Pine  Copper  Company  was  of  minimal  quality  and  had 
AC  voltages  as  large  as  50  volts  induced  on  it.  In  addition, 
the  transmission  of  higher  frequencies  was  severely  impaired. 

This  difficulty  was  overcome  by  installing  a new  data  line 
of  Spiral-4,  provided  by  the  Albuquerque  Center,  from  the 
subsurface  vaults  to  the  surface  recording  facility. 

Surface  installation.  During  this  report  period,  the  follow- 
ing alterations  of  the  White  Pine  surface  installation  have 
been  effected:  (1)  installation  of  the  uninterruptible  power 
system,  (2)  transferral  to  the  TG-120  timing  system  and  RV 
301  helicorders,  (3)  initiating  operation  of  the  Astrodata 
data  logger. 

The  Astrodata  data  logger  has  been  put  into  operation  with 
the  assistance  of  Albuquerque  personnel.  Several  integrated 
circuit  cards  were  faulty  and  had  to  be  replaced. 

Several  unavoidable  problems  arose/ all  of  which  were 

i 

eventually  surmounted  but  which  are  mentioned  in  the  follow-  I 

ing  extract  from  the  Technical  Report  1 April  - 30  September  i 

1975.  Station  installers  should  never  be  discouraged. 


Difficulties  encountered.  A major  mine  collapse  occurred 
on  August  11,  1975,  one  thousand  feet  to  the  southeast  of 
the  HGLP  vaults.  The  initial  collapse  was  300  feet  by 
300  feet  and  ran  to  the  surface.  At  approximately  the 
same  time,  the  east-west  instrument  began  having  signifi- 
cantly large  transients  which  have  decreased  in  number 
and  amplitude  with  time  but  still  persist.  The  north- 
south  instrument  has  displayed  the  same  characteristic 
but  to  a much  lesser  degree.  Since  the  Initial  collapse, 
the  cave-in  has  increased  in  size  to  300  feet  by  600 
feet  so  the  instability  persists.  Whether  the  instability 
of  the  east-west  instrument  is  due  to  the  collapse  or  an 
actual  instrumental  problem  is  not  known  at  this  time. 

If  the  instability  continues  when  the  vault  area  returns 
to  ambient  conditions,  the  east-west  instrument  will  have 
to  be  replaced.  The  collapse  poses  no  direct  threat  to 
the  HGLP  vaults. 

On  August  14,  1975,  a crane  operator  drove  through  the 
new  Spiial-4  data  line.  A splice  in  the  line  had  to  be 
effected  at  a pole  midway  between  the  vaults  and  record- 
ing area. 

Finally,  on  August  20,  1975,  the  WPM  observatory  was  shut 
down  by  high  voltage  from  a strike  or  near  strike  by 
lightning.  The  extent  of  the  damage  is  not  yet  completely 
rectified.  Three  discriminators  suffered  extensive  damage 
and  are  currently  being  repaired  by  Geotech.  The  mixer 
card  for  the  VCO's,  was  also  severely  damaged  and  had  to 
be  replaced. 

As  of  the  fall  of  1975,  the  station  was  recording  the 
velocity  signals  on  helicorders  at  a gain  of  approximately 
16K  and  at  full  gain  on  an  Astrodata  data  logger. 

During  the  period  October,  1975  through  March,  1976, 
the  station  was  intermittently  operational.  The  vertical 
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was  operational  only  about  65%  of  the  time  in  the  analog  mode 
and  the  horizontal  about  85%.  Much  of  the  downtime  was  due 
to  the  lightning  strike  discussed  earlier.  The  Astrodata 
data  logger  was  down  more  frequently.  Considerable  trouble 
was  encountered  with  large  random  transients  on  the  hori- 
zontals. 

In  the  Slimmer  of  197  6,  the  station  was  recalibrated 
and  put  into  routine  operation.  With  assistance  from  the 
Albuquerque  Seismological  Center,  difficulties  with  the  Astro- 
data  recorder  and  backup  power  systems  were  corrected.  In 
addition,  the  source  of  the  transients  on  the  horizontals, 
which  appeared  with  the  installation  of  the  Seismological 
Center  instruments,  was  located  and  corrected.  Insulating 
breaker  bars  of  pressed  phenolic  fiber  had  absorbed  enough 
moisture  to  cause  a periodic  short  across  one  or  more  of  the 
coils  of  both  horizontal  instruments.  The  breaker  blocks 
were  dried  and  desiccant  added  to  the  tanks  before  resealing. 

A series  of  calibration  response  curves  were  also  run  until 
the  response  curves  for  all  three  components  matched  the 
designated  HGLP  response  curve  to  within  5%  for  the  period 
range  10  to  300  seconds.  It  was  found  necessary  to  replace 
the  Ithaco  amplifier  to  match  the  response  curve  for  the 
vertical  instrument.  At  the  end  of  July,  all  instruments 
were  in  full  operation.  The  gains  on  the  helicorders  were 
set  at  30K. 

In  October  of  1976,  it  was  decided  to  close  the  station 
pending  decisions  as  to  the  feasibility  of  telemetering 
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the  signals  to  Ann  Arbor.  It  was  clear  at  this  time  that 
the  station  was  in  an  excellent  location.  Long  period  noise 
was  low  and  the  geographic  location  with  respect  to  inter- 
esting wave  paths  was  excellent;  however,  the  remoteness  of 
the  location  made  servicing  difficult  and  it  was  necessary 
to  depend  on  local  help  for  record  and  tape  changing. 

The  surface  recording  facility  was  removed  at  this  time 
and  moved  to  Ann  Arbor.  The  seismometers  and  associated 
equipment  underground  were  left  until  the  spring  of  1978 
when  it  was  clear  that  telemetry  was  not  feasible  in  the 
immediate  future. 

The  White  Pine  area  is  clearly  a prime  location  for  a 
high-quality  station  but  routine  operation  will  be  possible 
only  with  a station  engineered  for  minimum  maintenance  and 
with  the  data  telemetered. 

A considerable  amount  of  useful  information  regarding 
HGLP  installation  problems  and  their  solutions  was  obtained 
in  the  course  of  the  project.  Much  of  this  material  is  to 
be  found  in  the  technical  reports  and  has  been  shared  with 
the  Albuquerque  Seismic  Center.  It  is  hoped  that  eventually 
there  will  be  a high-quality  long  period  station  in  the 
Great  Lakes  region. 
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ABSTRACT 


Determinations  of  seismic  moment  have  been  made  for 
12  events  occurrinq  in  1974-1975  using  an  expanded  version  of 
the  AR  parameter  method  applicable  to  the  WWSSN  instruments. 

This  method  sums  the  energy  in  the  entire  surface  wave  train 
from  events  of  magnitude  4.0<M  <8.0,  recorded  over  predominantly 
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oceanic  paths  of  up  to  17,000  km  epicentral  distance  on  WWSSN 

long  period  three-component  instruments.  The  AR  values  are 

calibrated  directly  to  M.  The  moment  values  calculated 

include,  the  Utah-Idaho  March  28,  1975  event,  + 6.0,  with 

25 

new  moment  1.2  x 10  dyne-cm;  the  Leeward  Is.  October  8,  1974 

27 

event,  = 7.5,  with  moment  1.0  x 10  dyne-cm;  the  Peru 
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October  3,  197  4 event,  = 7 . 8,  with  moment  1.2  x 10  dyne-cm; 

the  Panama  July  13,  1974  event,  M = 7.3,  with  moment  6.0  x 
2 6 

10  dyne-cm;  the  Kermadec  July  2,  1974  event,  = 7.2, 

27 

with  moment  1.6  x 10  dyne-cm;  and  Solomon  Is.  February  1, 

27 

1974  event,  = 7.1,  with  moment  1.4  x 10  dyne-cm.  The 
method  is  shown  to  give  a factor  of  two  error  for  continental 
crustal  events,  and  a factor  of  three  error  for  oceanic 
crustal  events.  The  method  typically  requires  a much  smaller 
set  of  stations  than  in  spectral  amplitude  or  time  domain 
fitting  moment  estimates.  Only  an  approximate  knowledge  of 
the  recording  instrument  response,  path  type,  and  depth  are 
needed  for  application  of  the  method. 


INTRODUCTION 


Since  the  wide  acceptance  of  the  seismic  moment  (M^) 
as  a way  to  characterize  the  strength  of  an  earthquake 
source,  many  studies  have  been  conducted  to  calculate 
values  for  events  all  over  the  world  (Tsai  and  Aki,  1969; 
Kanamori,  1971;  Hanks  and  Wyss,  1972;  Mikumo,  1973; 

Langston,  1976;  Burdick  and  Mailman,  1976) . These  studies 
have  employed  a variety  of  methods  but  most  have  in  common 
the  need  for  digitized  seismograms  from  many  stations  for 
each  event  and  expensive  computer  fitting  routines  to 
determine  the  radiation  pattern. 

The  AR  parameter  method  was  introduced  by  Brune,  Espinosa, 
and  Oliver  (1963) , and  they  defined  the  AR  parameter  as  the 
sum  of  the  area  under  the  envelope  of  the  surface  wave 
amplitudes  recorded  on  long  period,  three-component  Press- 
Ewing  instruments.  In  Wyss  and  Brune  (1968,  1971)  the  AR 
value  was  related  directly  to  spectral  density, 

00 

AR  Lim  f f (t)  e dt 

Oj-vQ  -00 

They  assumed  that  the  is  proportional  to  the  spectral 

level  at  low  frequencies  as  predicted  by  seismic  source 

theory,  which  in  turn  relates  the  to  the  AR  parameter. 

One  Press-Ewing  long  period  instrument  was  calibrated  by 

measuring  AR  for  a set  of  events  in  the  Southern  California 

region  that  had  values  determined  by  other  methods.  Then 

the  calibration  curve  was  used  to  calculate  M for  other 

o 

events  in  that  same  locality.  These  workers  found  a direct 
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proportionality  of  and  the  AR  value  for  a set  of  small 
events  recorded  over  short  epicentral  distances.  They 
claimed  that  accurate  moments  could  be  found  without  applying 
radiation  pattern  corrections  since  adding  the  Love  and 
Rayleigh  surface  wave  trains  would  smooth  the  radiation  pattern 
effect. 

Several  recent  investigations  have  also  shown  the  method 
to  give  sufficiently  accurate  results  for  large  collections 
of  small  to  medium  magnitude  events  when  calibrated  to  a 
given  instrument  and  events  in  the  San  Andreas  transform 
fault  zone  (Hanks,  Hileman,  and  Thatcher,  1975)  and  in  the 
New  Zealand  Rift  zone  (Gibowicz  and  Hatherton,  1975) . 

In  light  of  the  potential  of  the  AR  parameter  method  to 
calculate  more  easily  and  routinely  than  other  methods, 
this  paper  develops  a single  AR  parameter  versus  calibration 
curve  for  20  long  period  WWSSN  stations  of  similar  tectonic 
settings,  recording  surface  waves  over  predominantly  oceanic 
paths,  epicentral  distances  up  to  17,u00  km  and  magnitudes 

4.0<M  <8.0.  The  calibration  curve  is  derived  from  a well 

s 

studied  set  of  events  with  known  and  radiation  patterns 
and  is  applied  to  12  unknown  events  occurring  in  1974  and 
1975.  The  effects  of  radiation  pattern  on  the  AR  value 
are  investigated  in  detail  using  the  Love  and  Rayleigh  wave 
spectral  density  formulations  of  Harkrider  (1970) . A 
theoretical  AR  vs  curve  is  derived  using  the  average 
spectral  theory  of  Geller  (1976)  . 
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DATA 


The  calibration  earthquakes  are  listed  in  Table  1 and 
are  from  a compilation  of  sources  by  Kanamori  and  Anderson 
(1975)  and  from  calculations  of  Gibowicz  and  Hatherton  (1975) . 
The  table  contains  several  moment  estimates  for  some  of  the 
events  aad  the  method  used.  The  unknown  data  set  for 
earthquakes  occurring  in  1974-1975  is  in  Table  2. 

Stations  were  chosen  for  measurement  that  had  three 
components  of  the  long  period  WWSSN  instruments  recorded, 
with  Tq  = 15  sec  and  T^  = 100  sec.  Each  event  has  AR  values 
from  several  stations;  the  widest  range  of  epicentral  distances, 
from  422  to  17,000  km,  and  widest  azimuthal  coverage  possible. 
These  st^  ons  are  listed  in  Table  4 for  each  event.  The  paths 
to  each  station  are  restricted  to  predominantly  ocean  paths 
to  simplify  attenuation  corrections  later.  Gains  range  from 
3000  to  375  and  paper  speeds  from  15  mm/min  and  30  mm/min. 

The  AR  values  are  normalized  to  a gain  of  750  and  paper  speed 
of  30  mm/min  on  a standard  size  WWSSN  seismogram. 

The  envelope  of  the  surface  waves  was  drawn  on  Itek 

2 

seismogram  reproductions,  and  converted  to  AR  in  cm  . This 
will  be  called  AR  (A) . The  larger  than  original  size 
seismographs  facilitated  the  separation  of  successive  traces 
on  the  record,  which  could  be  a significant  source  of  error 
in  the  measurement  technique  for  the  large  earthquakes  in 
this  study.  With  care,  even  for  8.0,  the  surface  wave 

envelope  may  be  measured  with  errors  of  less  than  10%. 
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compared  to  the  radiation  pattern  variations. 

Figure  1 shows  several  AR  envelopes  drawn  in  this  study 
and  some  of  the  techniques  used  for  handling  special  problems. 
These  procedures,  some  derived  from  Brune  et  a]^.  (1963)  , are 
described  in  detail  in  Williams  (1977),  and  tended  to  reduce 
errors  and  time. 

The  effect  of  aftershocks  or  events  occurring  in  the 
wave  train  of  the  main  event  is  a major  problem  when  measuring 
the  AR  value  for  larger  events.  If  there  are  several 
significant  aftershcoks  within  2 magnitude  units  of  the  main 
shock,  an  attempt  is  made  to  identify  by  inspection  the 
arrival  of  these  events  within  the  wave  train,  measure  the 
area  under  the  envelope  to  the  20  s wave  arrival  time,  and 

subtract  this  area  to  form  a minimum  and  maximum  AR  value. 

2 

The  resulting  AR(A)  values  in  cm  for  the  calibration 
and  unknown  data  sets  are  plotted  vs  epicentral  distance  in 
Fig.  2.  These  will  serve  as  the  basis  for  deriving 
corrections  for  AR(A)  with  distance.  ..e  will  show  that  it  is 
not  necessary  to  restrict  the  measurements  to  one  station 
and  calibrate  it  separately  to  find  a station  correction 
as  Brune  £t  (1963)  , Wyss  and  Brune  (1968)  , Gibowicz  and 

Hatherton  (1975)  and  Hanks,  Hileman,  and  Thatcher  (1975)  have 
done.  Only  one  station  gave  consistently  high  AR  values, 
the  South  Pole  station  SPA.  This  factor  of  about  1.7 
remained  fairly  constant  along  the  magnitude  range.  In 
Fig.  2,  the  station  SPA  data  has  been  reduced  to  and  plotted 
with  the  other  stations. 
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Change  of  AR(A)  with  Epicentral  Distance 


The  success  of  the  AR  method  depends  on  how  accurate  a 
correction  of  AR(A)  with  epicentral  distance  one  can  make. 

Brune  et  al.  (1963)  simplified  this  procedure  by  making  the 
following  assumption  at  each  distance:  one  can  characterize 
the  decay  or  growth  of  the  AR  envelope  by  one  period,  and 
as  the  AR  is  measured  at  larger  distances  this  representative 
period  increases  because  the  higher  frequencies  become 
attenuated. 

In  Fig.  3,  Brune ' s five  step  model  is  shown  for  the 
representative  period  at  various  distances.  Since  we  have  a 
wide  range  of  magnitude  curves  A,B,  and  C have  been  adopted 
as  the  models  for  small,  medium,  and  large  events,  assuming 
in  general  that  larger  magnitude  events  will  show  longer  period 
contributions  to  the  surface  wave  train  than  the  next  smaller 
group.  Curve  A is  below  Brune 's  model  because  the  Press-Ewing 
long  period  instruments  (Tq  = 30  s,  = 100  s)  he  used  have 
a slightly  larger  period  response  peak  than  the  present 
setting  of  the  long  period  WWSSN  instruments  (T^  = 15  s, 

Tg  = 100  s) . Several  models  have  been  tried  in  this  analysis 
and  those  that  best  fit  the  measured  AR(A)  data  in  Fig.  3 
are  shown  as  models  A,  B,  C. 

Four  main  corrections  will  be  applied  to  the  AR(A)  that 
depend  on  period  or  epicentral  distance;  geometric  spreading, 
instrument  response,  attenuation,  and  spreading  due  to 
dispersion  of  the  wave  train.  No  corrections  for  crustal 
response  or  source  spectral  shape  will  be  applied  to  the  AR(A) 
at  this  time. 
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The  geometric  spreading  used  is  that  adopted  by  Brune 
et  al.  (1963),  a factor  of  l//sinA°  applied  to  the  AR{A) . The 
instrument  response  is  based  on  the  magnification  curve  for  a 
WWSSN  long  period  instrument  at  750  gain.  The  attenuation 
corrections  for  the  Pacific  Ocean  paths  are  based  on  Mitchell 
et  al.  (1976) . These  coefficients  are  used  to  correct  for 
Atlantic  and  Indian  Ocean  paths  as  well.  The  Love  and 
Rayleigh  coefficients  are  averaged  for  each  period.  The  coef- 
ficients are  applied  via  the  relation  exp- (y^^  + Yj^/2)A,  where 
and  are  the  Love  and  Rayleigh  attenuation  coefficients 
for  the  period  predominant  at  the  distance  A of  the  station, 
according  to  the  period  model  A,  B,  or  C. 

After  Brune,  the  change  in  the  AR(A)  due  to  spreading 

1/2 

of  the  train  is  a factor  of  T where  T is  the  duration  of 
the  wave  train.  Brune  assumed  that  for  his  set  of  small 
earthquakes,  the  range  of  group  velocities  was  constant, 
therefore  the  duration  grew  linearly  with  distance.  But 
in  this  study,  oceanic  paths  are  used  and  the  range  of  group 
velocities  will  differ  between  curves  A,  B,  and  C.  The  upper 
limit  of  the  goup  velocity  is  set  at  4.4  km/sec  which  is  the 
Love  wave  oceanic  group  velocity  for  the  periods  of  interest 
in  this  study.  The  group  velocity  lower  limit,  though 

will  change  rapidly  with  period  due  to  the  steep  Rayleigh  wave 
oceanic  dispersion  curve  and  the  dispersion  factor  will 
depend  on  how  fast  the  lower  velocities  become  neglected. 

This  factor  is  D(A)  = [A  d/Unjin  “ 1/4.4)]^/^  and  Fig.  4 
shows  the  change  in  duration,  D(A)  of  the  wave  train  due  to 
the  effect  of  group  velocity  dispersion,  versus  distance,  A. 
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The  curve  corresponding  to  the  A model  shows  a steeper 
slope  than  Brune's  as  the  velocity  range  is  increased  through 
the  Airy  phase,  then  flattens  as  the  low  group  velocities 
become  neglected.  B'  shows  an  intermediate  model  where  the 
Airy  phase  velocities  are  neglected  at  about  7000  km.  The 
curve  corresponding  to  model  B neglects  these  low  velocities 
more  gradually  than  B',  starting  at  about  2000  km.  The  C 
curve  has  a nearly  constant  velocity  range. 

If  we  generate  decay  curves  of  AR(A)  with  distance  at 
this  point  and  correct  the  measured  data  sets  to  a standard 
distance,  the  small,  medium,  and  large  groups  of  events  will 
be  offset  from  each  other  because  the  starting  points  of  the 
three  models  A,  B,  and  C are  different  in  period  and  group 
velocity  range.  Several  standards  must  be  chosen  which  will 
make  the  AR(A)  corrected  to  a given  distance  consistent  between 
various  magnitudes  and  kinds  of  events. 

The  standards  are  as  follows: 

1.  All  AR(A)  will  be  corrected  to  the  value  at  3.8°, 
the  distance  at  which  dispersion  is  completely  developed 
for  oceanic  paths  (Gibowicz  and  Hatherton,  1975) . 

2.  Attenuation  and  instrument  response  characteristics 
will  be  normalized  to  those  at  20  s. 

3.  All  durations  of  dispersed  wave  trains  will  be 
scaled  to  the  standard  duration  from  a wave  train  with  group 
velocities  of  .5-4.5  km/sec. 

4.  A fourth  standard  is  implied  by  the  use  of  predominantly 
shallow  continental  crust  events  in  the  calibration  data  set. 


The  oceanic  crustal  response  is  generally  lower  than  the 
continental  crustal  response  for  a large  range  of  periods 
of  interest  in  this  study  and  a variety  of  shallow  source 
mechanisms  (Harkrider,  1970) . Thus,  for  events  occurring 
in  ocean  crust  settings,  another  scale  factor  calculated 
from  Harkrider 's  formulations,  is  applied  to  raise  the 
crustal  response  to  the  standard;  shallow  events  occurring 
in  continental  crust.  The  Gutenberg-Bullen  earth  model  is 
the  continental  crustal  model,  and  the  Harkrider  oceanic 
crust  model  is  chosen  to  represent  the  oceanic  tectonic 
settings . 

No  corrections  for  depth  are  made  in  this  study.  All 
earthquakes  are  assumed  to  be  40  km  or  shallower.  If  depth 
corrections  are  desired,  they  would  be  applied  to  the  combined 
standard  scale  factors. 

The  final  three  curves  based  on  representative  period 
models  A,  B,  and  C of  AR(A)  decay  with  epicentral  distance 
are  presented  in  Fig.  5 along  with  the  combined  standard 
scale  factors  for  each  model  and  crustal  type.  The 
theoretical  decay  curves  in  Fig.  5 show  three  different 
slopes  in  the  2000-10,000  km  range  which  can  be  used  to 
distinguish  the  three  groups  of  events;  the  A model  cxirve 
has  the  steepest  slope,  B an  intermediate,  and  C the 
smallest  slope.  Fig.  2 shows  that  these  slope  changes  do 
exist  in  the  data  by  comparing  the  group  of  numbered  earth- 
quakes 68/76,  68/146,  and  67/426  to  the  Turkey  '67  and  Iran 
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events.  As  predicted,  there  is  a definite  decrease  in 

slope  as  the  AR  values  become  larger.  The  first  break  point 

seems  to  come  at  M -v  6.0,  and  the  second  break  point  at  M 'v  7.0. 

= s 

This  will  serve  as  a general  guideline.  Events  with  < 6.0 

should  fit  closely  to  the  A model  decay  curve,  events  with  6.0  < 

M < 7.0  the  B model  decay  curve,  and  for  M >7.0  the  C model 
s s 

curve.  Regions  on  Fig.  2 outline  the  AR(A)  values  to  be 
applied  to  each  curve. 

There  are  some  exceptions  to  the  division.  The  Borrego 

Mountain  event  with  M =6.7  definitely  has  the  characteristics 

s 

of  the  A decay  curve.  But  this  event  has  an  exceptionally 
large  number  of  aftershocks  which  show  up  on  the  closer 
stations.  These  will  attenuate  quite  fast  as  the  epicentral 
distance  is  increased,  so  the  data  points  will  give  an 
apparently  steep  slope.  The  Borrego  Mountain  event  will  be 
fit  to  the  B curve  using  the  minimum  AR(A)  values. 

The  Solomon  Islands,  Leeward  Islands,  Kermadec  1,  and 
India  events  show  large  kinks  in  their  AR(A)  values.  These 
deviations  will  be  shown  in  the  next  section  to  be  due  to  the 
radiation  pattern  effect  for  these  events. 

Figures  6a,  b,  and  c show  the  observed  AR(A)  values  fit 
to  the  proper  theroetical  decay  curve.  The  earthquakes 
chosen  to  fit  each  decay  curve  are  labeled  on  the  plot. 


Azimuthal  Variation  of  the  AR.  _o  Parameter 

J « O 

To  Study  the  azimuthal  variation  of  the  AR-  -o  parameter, 

J • O 

plots  of  the  theoretical  surface  wave  radiation  pattern  for 
the  calibration  set  of  events  and  for  four  events  of  1974  will 
be  generated.  Focal  mechanism  solutions  have  been  determined 
by  previous  investigators  and  are  listed  in  Table  3.  The 
AR^  gO  values  are  drawn  on  the  plots  in  their  azimuthal  position 
to  gauge  the  extent  of  the  radiation  pattern  and  path  difference 
effects. 

Since  the  radiation  pattern  of  surface  waves  changes  with 
period,  several  periods  were  chosen  for  illustration  based 
upon  the  A,  B,  or  C representative  period  models.  For  the 
events  fit  to  the  A curve  in  the  reduction  of  the  AR{A) 
to  AR,  qO,  the  20  s radiation  pattern  is  compared  to  the  AR,  -o 

0*0  . j • O 

values,  for  the  B curve  events,  a 30  s pattern  is  chosen  and 
for  the  C curve  events,  a 45  s pattern. 

Since  the  AR,  „o  parameter  is  related  to  the  addition 
of  the  spectral  density  for  each  component,  it  will  be  compared 
to  the  spectral  far  field  solutions  for  Love  and  Rayleigh 
waves  from  Harkrider  (1970) . The  same  source  function 
used  by  Geller  (1976)  is  chosen,  M(w)/a),  in  which  he  develops 
average  theoretical  spectra  for  various  magnitudes  of  events. 

The  envelope  radiation  diagram  will  show  the  E^(a))  function, 

E^(co)  = R(a))  exp(Yj^-Yp)A  + |L(a));  (1) 
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where  iR(u))|  is  the  addition  of  the  horizontal  and  vertical 
Rayleigh  wave  responses  and  |L(a))  | is  the  Love  wave  response. 
These  values  are  calculated  from  Harkrider's  (1970)  formulations 
and  depend  on  crustal  model,  the  source  mechanism  and  spectrum. 
The  exponential  term  is  used  to  adjust  the  theoretical  envelope 
for  relative  differences  in  the  attenuation  of  Love  and 
Rayleigh  waves  over  oceanic  paths  which  will  affect  the  shape 
of  the  observed  AR  parameter  azimuthal  distribution. 

The  AR^  go  parameter  must  now  be  related  to  the  spectral 
densities.  The  steepest  descent  approximation  for  surface  waves 
relates  the  amplitudes  at  the  phase  period  in  the  frequency 
domain  to  the  group  period  in  the  time  domain  (Ben  Menahem, 

1961) . The  amplitudes  in  the  time  domain  are  modified  by  the 
slope  of  the  dispersion  curve  at  the  group  period  and  by 
the  phase  shifts  introduced  by  assuming  a finite  earthquake 
source  with  duration  time.  These  phase  shifts  have  been 
assumed  to  be  small  except  for  the  largest  events  on  record 
(Kanamori  and  Anderson,  1975;  and  Geller,  1976)  and  have 
been  neglected  in  this  study.  The  final  relation  is  from 
Williams  (1977)  , 


|r(u))  e"^R^  + L(to)  e~^L'^ 


—202 

A . M(a)).  n-10  “^^cm 


The  exponential  factor  is  the  attenuation  correction  adopted 
in  the  formulation  of  the  AR  decay  curves  with  epicentral 
distance,  an  average  of  the  Love  and  Rayleigh  coefficients. 
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The  scale  factor  corrects  for  the  standard  set  of 

parameters,  i.e.,  422  km  distance,  750  gain,  30  mm/min 
paper  speed,  attenuation  and  instrument  response  at  20  s. 
This  factor  also  scales  for  the  fact  that  the  AR-  „„  value 

^ • O 

is  calculated  from  a dispersed  wave  train.  Developing  the 
above  equation  further 

AR 

— — e (Yt  -Yr/2)A  = E-  (co)  sec/dyne  (2) 

M(a))  n ^ ^ 

Figure  7 gives  the  calculated  and  observed  envelope  func- 
tions from  equations  (1)  and  (2)  for  continental  crustal 
events  in  this  study  that  have  wide  azimuthal  coverage. 
Figure  8 shows  the  same  for  oceanic  crustal  events.  The 
observed  envelope  functions  measured  from  AR_  oo  values 
are  plotted  at  each  station's  azimuth. 

Stations  located  in  predicted  nodal  positions  rarely 
give  AR^  go  values  as  low  as  predicted,  but  Brune's  assump- 
tion that  adding  the  love  and  Rayleigh  waves  smooths  the 
radiation  pattern  effect  can  be  seen  to  work  nearly 
completely  only  for  strike-slip  events  in  continental 
crust.  The  Iran  and  Borrego  Mountain  events  show  little 
observed  or  predicted  azimuthal  variation  in  AR-  oo  because 
the  Love  and  Rayleigh  wave  response  is  nearly  equal.  Even 
the  India  earthquake,  more  nearly  a dip-slip  event,  shows 
only  a factor  of  2 variation.  Contrast  these  continental 
crustal  events  to  those  in  Figure  8,  where  the  dip-slip 
events  in  oceanic  crust  produce  nearly  a factor  of  3 
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variation  as  the  Love  wave  response  becomes  much  smaller 
than  the  Rayleigh  wave  response.  Dip-slip  events  concen- 
trate the  Rayleigh  wave  energy  near  90°  to  the  strike  of 
the  fault  plane.  These  azimuths,  45°  to  135°  from  the 
strike  of  the  fault  plane,  would  be  a favored  position  for 
a set  of  AR.,  oo  measurements  because  the  maximum  energy 

J • o 

from  both  strike-slip  and  dip-slip  events  would  be  recorded. 

Since  the  radiation  pattern  and  aftershock  activity 

have  such  a dominant  effect  on  the  AR(A)  data  fits  to  the 

theoretical  curves  in  Figs.  6a,  b,  and  c,  the  errors  are 

best  discussed  in  the  context  of  certain  classes  of  events 

rather  than  based  upon  the  measuring  procedures  for  each 

event.  These  procedures  contribute  less  than  10%  error  to 

the  AR(A)  value  for  all  magnitude  ranges. 

The  classes  of  events  which  have  the  best  fit  to  the 

AR(A)  data  to  theoretical  curve,  to  within  a factor  of  two, 

are  the  continental  crustal  events  and  small  oceanic  events 

of  M < 6.0,  both  without  aftershock  activity.  These  errors 
s 

are  caused  by  minor  path  variations,  the  presence  of 

multipathing  in  the  wave  train,  and  a small  radiation  pattern 

effect.  Continental  crustal  events  with  many  aftershocks, 

and  oceanic  events  with  M >6.0  and  little  aftershock 

s 

activity  fit  the  theoretical  curves  by  a factor  of  three. 

The  extra  source  of  error  for  the  continental  events  is  the 
contamination  of  the  wave  train  by  interfering  events,  and 
the  error  in  the  larger  oceanic  earthquakes  is  due  to  the 
increased  radiation  pattern  effect  for  all  sources  within 
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oceanic  crust.  A factor  of  four  error  is  assigned  to  the 
last  class  of  events,  the  large  oceanic  crustal  events 
with  much  aftershock,  activity.  The  list  of  values 
generated  by  different  techniques  for  the  calibration 
earthquakes  in  Table  1 shows  that  this  is  an  acceptable 
kind  of  error  for  these  calculations. 

Calibration  Curves 

All  the  data  was  reduced  to  the  standard  distance  to 
derive  the  AR^  go  value  using  the  appropriate  decay  curve 
and  standard  corrections.  These  values  are  plotted  for 
the  calibration  earthquakes  in  Fig.  9 versus  the  accepted 
for  these  events.  Error  bars  have  been  drawn  for  the 
AR^  go  and  values.  In  general,  the  AR^  go  parameter 
agrees  quite  well  within  the  error  bars  and  forms  a 
consistent  calibration  curve  over  a wide  range  of  moment 
values . 

The  relations  for  the  calibration  curve  are 

Log  = Log  AR^  goCm^  + 22.0  AR^  go  < 3 x 10^  cm^ 

Log  = 1.55  Log  AR^  goCm^  +19.5  AR^  go  > 3 x 10^  cm^ 

The  ARg  go  data  from  the  unknown  set  of  events  in  1974-1975 
are  converted  to  using  the  above  relationships.  To  arrive 
at  a final  estimate,  the  stations  are  averaged.  This  data 
is  in  Table  4. 

The  strategic  location  of  station  SPA  plus  its  low 
gain  setting  (375)  makes  it  especially  suited  to  record 


large  events  in  the  Pacific.  This  station  was  the  only  one 
to  give  consistently  higher  AR(A)  measurements.  Now  part 
of  the  calibration  and  unknown  events  can  serve  as  a cali- 
bration for  this  station.  This  data  is  shown  in  Fig.  10. 

The  tectonic  setting  of  station  SPA  is  quite  different  from 
that  of  the  rest  of  the  stations;  SPA  sits  well  into  the 
continent  of  Antarctica  while  most  of  the  other  stations  sit 
near  the  oceanic-continental  edge.  These  variations  in 
tectonic  setting  do  not  give  a correction  of  more  than  a 
factor  of  two  for  SPA  and  for  most  of  the  magnitude  range 
the  factor  is  less.  In  light  of  this,  the  assumption  that 
stations  with  similar  paths  from  the  epicenter  could  be 
used  without  calibrating  each  separately  seems  reasonable. 

The  relations  for  the  SPA  calibration  curve  are 

Log  = Log  AR^  go  + 21.8  AR^  go  £ 3.5  x 10^  cm^ 

Log  = 1.63  Log  AR^  go  + 18.9  ARg  go  > 3.5  x 10**  cm^ 

Theoretical  and  Observed  AR-  -o  vs  M Relations 

.3*0  O 

A theoretical  AR,  -o  vs  M calibration  curve  can  be 

J • o O 

generated  for  various  frequencies  from  relation  (2).  Using 
the  average  M(a))  defined  by  Geller  for  a stress  drop  of 
50  bars,  and  an  average  value  for  the  20  s Rayleigh-Love 
envelope  response  function,  the  theoretical  AR,  ^o  vs  M 
curve  with  constant  crustal  response  is  presented  in  Fig.  11. 
This  set  of  curves  will  move  to  the  right  for  higher  stress 
drop  than  50  bars,  and  to  the  left  for  lower  stress  drop. 
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( 


These  curves  reflect  the  spectral  shape  change  with 
increasing  from  source  theory.  Fig.  12  shows  the 
addition  of  variable  crustal  response.  The  two  effects 
of  crustal  response  with  period  and  spectra  shape  change 
with  increasing  were  the  only  two  corrections  left  out 
of  the  AR^  go  snd  so  must  be  applied  to  the  theoretical 
curves . 

Fig.  13  compares  the  observed  calibration  curve  and 

the  envelope  of  selected  theoretical  curves.  This  envelope 

increases  in  period  with  increasing  which  is  the 

assumption  adopted  in  drawing  the  predominant  period  models 

A,  B,  and  C.  The  fit  of  the  two  curves  is  quite  good,  and 

26 

the  steeper  slope  for  above  3 x 10  dyne-cm  is  inter- 
preted as  a source  and  crustal  response  effect.  This 
comparison  establishes  the  calibration  curve  in  the  AR^  go 
parameter  method  to  a true  static  in  the  context  of 
seismic  source  theory. 

CONCLUSION 

A single  calibration  curve  for  relating  the  AR 

parameter  to  has  been  found  which  is  applicable  to  20 

long  period  WWSSN  stations  recording  surface  waves  from 

shallow  earthquakes  over  predominantly  oceanic  paths, 

epicentral  distances  up  to  17,000  km  and  magnitudes  up  to 

M =8.0.  The  relation  is 
s 

Log  = Log  AR^  goCm^  + 22.0  AR^  go  1 3 X lO^cm^ 

Log  = 1.55  Log  AR^  goCm^  + 19.5  AR^  go  > 3 x lO^cm^ 
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This  curve  relates  well  to  that  predicted  by  seismic 


source  theory  where  is  interpreted  as  a true  static 

M value, 
o 

The  procedure  is  quick  and  easy,  requiring  only  an 
approximate  knowledge  about  the  recording  instrument 
response,  path  type,  and  depth.  For  continental  events, 
few  stations  may  be  needed  to  calculate  a within  a 

I 

factor  of  2 error,  for  the  AR  parameter  shows  little 

radiation  pattern  variation.  For  oceanic  events,  the 

radiation  pattern  effect  on  the  AR  parameter  becomes  more 

pronounced,  and  a set  of  stations  azimuthally  distributed 

are  needed  to  calculate  M values  with  factor  of  3 error. 

o 

Typically,  a much  smaller  set  of  stations  is  required 
for  the  AR  parameter  than  in  spectral  amplitude  or  time 
domain  estimates  of  moment  because  the  AR  value,  being  a 
sum  of  many  periods,  is  less  sensitive  to  path,  station 
and  source  differences.  Thus,  can  be  calculated  more 
routinely  than  with  the  methods  predominantly  in  use. 


This  research  was  supported  by  the  Advanced  Research 
Projects  Agency  of  the  Department  of  Defense  and  was 
monitored  by  the  Air  Force  Office  of  Scientific 
Research  under  Contract  No.  F44620-73-C-0060. 
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Event  Type;  A - small  oceanic  crust,  B - continental  crust,  C - large  ocean  crust 
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Events  with  Unknown 

o 


Event 

Date 

Depth 

(Km) 

M  *  ** 
8 

Solomon  Island 

February  1,  1974 

12 

6.2 

7.1 

Szechwan  Province, 
China 

May  10,  1974 

17 

00 

• 

in 

6.8 

Mariana  Island 

May  11,  1974 

26 

6.2 

5.9 

Molucca  Pass. 

May  11,  1974 

43 

6.1 

6.0 

South  India  Ocean 

June  25,  1974 

20 

6.1 

6.7 

Kermadec  1 

July  2,  1974 

33 
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77 

6.0 

6.6 

Panama 

July  13,  1974 

12 
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7.3 

Kurile  Island 
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5 
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6.7 
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October  3,  1974 

9 
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7.8'^ 
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30'"'^ 
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7.4 

Utah- Idaho 

March  28,  1975 

5 
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* ISC  Bulletin 

**  NEIS 

+ Stewart  (1978) 

++  TonUDlin  and  Aspinall  (1975) 
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Table  3 


Radiation  Pattern  Input  Data 


Strike 


Event  deg. 

India  16 

Borrego  Mt.  315 

Iran  275 

Tokachi-Oki  156 

Kermadec  1 184 

Panama  245 
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45 
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Langston  (1976) 

Burdick  & Mellman  (1976) 
Niazi  (1969) 

Kanamori  (1971) 

Williams  (1977) 

Williams  (1977) 

Williams  (1977) 


Solomon  Is 


Table  4 


Event 

Curve 

A,B,C 

!!o 

sta- 
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X lO^cm^ 

AR3  go  Method 
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0 

X lO^e 

Average 

Mo 
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Radiatior 
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25 
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27 
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AKU 

91.4 

17 
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C 
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1.1x10^® 

2i 
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C 
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2f 
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C 
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^ _ 
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p. 
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applied. 
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Figure  5. 
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Figure  6c  Fit  of  large  events  to  C model  decay  curve  from 
Figure  5.  Note  increase  in  spread  of  data  due 
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ABSTRACT 


A new  model  to  predict  Rayleigh  wave  group  velocity  dispersion  (per- 
iods 30  through  100  seconds)  through  the  ocean  basins  is  presented. 

It  is  based  on  a 5°x  5°  integration  of  weighted  Rayleigh  group  vel- 
ocities for  the  2592  geographic  grid  elements  constituting  the  earth's 
surface.  Percentages  (to  the  nearest  10%)  of  element  coverage  by 
each  of  twenty  tectonic  divisions  provide  the  weighting  factors.  Ray- 
leigh group  velocities  specific  to  each  of  the  twelve  continental 
and  one  oceanic  geological  division  were  compiled  from  the  literature. 
'Rayleigh  velocities  for  the  remaining  seven  oceanic  divisions  were 
determined  by  partitioning  empirically-derived  functional  relationships 
between  sea  floor  age  and  Rayleigh  group  velocities. 

Dispersion  analyses  for  seventeen  random  paths  through  the  Pacific 
Ocean,  four  for  the  Indian  Ocean  and  six  for  the  Atlantic  Oceans  are 
compared  with  the  model -predicted  dispersion  curves.  Observed-minus- 
predicted  velocity  residuals  for  the  Pacific  are  significantly  less 
than  the  uncertainty  limits.  Systematic  positive  residuals  occur  in 
basins  of  the  other  major  oceans,  however.  The  largest  anomaly  (0.4 
km/sec)  occurs  in  the  Brazil  Basin.  The  failure  of  the  Pacific-based 
model  to  predict  dispersion  accurately  in  these  other  basins  suggests 
that  there  may  be  fundamental  differences  in  crust  and  upper  mantle 
structure  among  the  world's  oceans. 

Since  petrological  differences  among  the  oceans  are  apparently 
negligible,  perturbations  of  'typical'  oceanic  geotherms  by  lower 
thermal  influx  (possibly  reflected  in  lower  spreading  rates)  and/or 


thermal  refractions  into  colder  adjacent  continental  regions  may 
account  for  the  observed  velocity  differences.  Lower  temperatures 
would  cause  a reduction  of  partial  melting  which  would  result  in 
faster  velocities  through  the  region. 
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Introduction 

Analysis  of  earthquake  and  explosion  generated  Love  and  Rayleigh 
mode  surface  waves  has  contributed  significantly  to  the  elucidation 
of  the  crust  and  upper  mantle  structure  of  the  earth.  Although  the 
dispersion  of  surface  waves  is  determined  by  the  average  properties 
of  the  crust  and  mantle  over  a large  depth  range,  sampling  phase  and 
group  velocities  of  dispersed  wave  trains  at  a variety  of  periods  has 
revealed  fundamental  structural  differences  between  continents  and 
oceans  (Oliver,  1962).  In  addition,  variations  of  observed  continen- 
tal surface  wave  dispersion  character  have  been  correlated  convincing- 
ly with  geological  provinces  (North,  1976;  Goncz  and  Cleary,  1976,  Fix, 
1975;  Teng,  1974;  Biswas  and  Knopoff,  1974;  Adams,  1974;  Wu,  1972; 

James,  1971;  Bolt  and  Niazi,  1964;  Brune  and  Dorman,  1963;  Brilliant 
and  Ewing,  1954).  Similar  variations  of  oceanic  surface  wave  disper- 
sion character  were  recognized  before  sea  floor  spreading  was  introduced 
(Tarr,  1969;  Sato  and  Santo,  1969;  Santo,  1966a,  1966b,  1962,  1961a, 
1961b,  1960;  Papazachos,  1964;  Brune  and  Major,  1962;  Kovach,  1961; 
Shurbet,  1960;  Berckhemer,  1956;  Oliver,  Ewing  and  Press,  1955;  Jardetsky 
and  Press,  1953;  Wilson  and  Baykal,  1948),  but  without  the  plate  tec- 
tonics and  ocean  floor  evolution  concepts,  a rational  differentiation 
of  observed  dispersion  into  provinces  was  elusive. 

The  fundamental  discovery  that  the  sea  floor  is  created  at  ocean 
ridges  and  evolves  both  spatially  and  temporally  provided  the  key  to  a 
rational  catagorization  of  the  diverse  dispersion  characteristics  pre- 
viously observed  in  the  ocean  basins.  The  development  of  the  ocean  floor 
spreading  framework  from  marine  magnetic  anomalies  and  deep  sea  drilling 
has  enabled  a chronologically  detailed  mapping  of  the  world's  oceans 
(Pitman,  Larson,  and  Herron,  1974;  Herron,  Pitman,  and  Dewey,  1974; 
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Cooper,  Marlow  and  Scholl,  1976;  Johnson,  Powell  and  Veevers,  1976). 
Equipped  with  this  new  information  on  the  age  of  the  ocean  floor,  and 
with  improved  techniques  of  dispersion  analysis  (Knopoff,  1972;  Landis- 
man,  Dziewonski  and  Sato,  1969),  Forsyth  (1976,  1975),  Yoshii  (1975), 

Leeds  (1975),  Leeds  Knopoff  and  Kausel  (1974a, b)  and  Knopoff,  Schlue  and 
Schwab  (1970)  have  demonstrated  the  correlation  between  ocean  floor  age 
and  variations  in  both  Love  and  Rayleigh  dispersion  in  the  Pacific  basins. 

The  research  discussed  in  this  paper  is  an  extension  of  the  character- 
ization of  Rayleigh  group  velocity  dispersion  through  ocean  basins.  An 
empirical  model  based  on  sea  floor  age  used  to  predict  Rayleigh  group 
velocity  dispersion  for  any  path  through  ocean  basins  is  presented  and 
discussed.  The  principal  model  assumption  is  that  Rayleigh  group  vel- 
ocities are  definable  and  constitute  well-behaved  functions  in  the  tectonic 
evolutionary  context  of  sea  floor  spreading.  Global  characterization  of 
Rayleigh  group  velocities  in  a 5°x5°  tectonic  framework  requires:  first, 
a regionalized  encoding  of  global  tectonic  units  (e.g.  Archean  shield, 

0-5  million  year  old  oceanic  crust,  etc.),  second,  an  assignment  of  a 
characteristic  dispersion  through  each  tect'^nic  unit,  and  third,  a merging 
of  the  5°x5°  encoded  tectonics  with  the  tectonic-specific  group  velocities 
to  obtain  a global  representation  of  the  period-specific  Rayleigh  group 
velocity  fields.  Inversion  of  the  predictive  velocity  model  to  obtain  a 
specific  crust  and  upper  mantle  velocity  structure  is  not  intended. 

CONSTRUCTION  OF  THE  GLOBAL  RAYLEIGH  GROUP  VELOCITY  MATRIX 

A fundatnental  tectonic  data  matrix  was  created  by  dividing  the 
surface  of  the  earth  into  a grid  of  2592  5°x5°  elements  and  encoding 
the  tectonic  information  for  each  element  derived  from  maps  and  pub- 
lished research.  Although  this  grid  does  not  yield  elements  of  equi- 
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valent  area  (equatorial  -v-  308,715  km^,  polar  -^13,478  km^),  it  does 
provide  a convenient  system  for  subsequent  computational  manipulation. 

Twelve  continental  and  eight  oceanic  divisions  were  used  for  the 
matrix  encoding  (See  Table  1.).  Although  many  of  the  tectonic  divisions 
might  be  subdivided  further  for  specific  regions,  the  selected  larger 
scale  designations  could  be  universally  recognized  and  were  deemed  suf- 
ficient for  this  global  study. 

The  sources  of  the  tectonic  data  are  listed  in  the  appendix.  Where 

possible,  detailed  tectonic  data  are  used  to  supplement  or  to  supersede 

older  and/or  more  generalized  data.  The  Indian  Ocean  data  of  Johnson, 

Powell  and  Veevers  (1976),  superseded  the  data  of  McKenzie  and  Sclater 

(1971),  for  example.  This  was  the  largest  scale  revision.  More  fre- 
quently, tectonic  revisions  involved  only  one  or  two  element  adjustments. 

The  tectonic  data  matrix  is  amendable  as  new  data  become  available.  The 

latest  update  used  in  this  study  was  March,  1976. 

Generally,  continental  tectonic  data  are  available  in  abundance.  In- 
accessabil ity  (central  Asia),  or  ice  cover  (Greenland  and  Antarctica) 
may  preclude  detailed  analysis  of  some  regions,  but  a reasonable  ap- 
proximation of  the  tectonics  could  still  be  effected.  Significant  gaps 
in  the  coverage  of  the  ocean  floor  do  exist,  however.  Detailed  tectonic 
data  for  the  Arctic  Ocean  Basin,  for  example,  are  sparse.  In  this  case, 
a postulated  model  for  Arctic  evolution  proposed  by  Herron,  Dewey  and 
Pitman  (1974)  was  used  to  fill  the  gap.  Similarly,  gaps  in  the  ocean 
floor  coverage  in  the  equatorial  Atlantic  and  Atlantic-Indian  Basin  are 
problematic.  In  these  instances,  where  no  direct  data  exist,  tectonic 
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patterns  were  inferred  on  the  basis  of  immediately  adjacent  element  cov- 
erage and  other  tectonic  information  such  as  seismicity  and  transform 
fault  offsets.  Finally,  the  marginal  basins  often  are  poorly  dated.  The 
majority  of  them  .however,  have  deep  sea  drilling  logs  and/or  subsequent 
detailed  studies  to  reduce  the  uncertainty  of  tectonic  assignment.  These 

V 

ill -defined  ocean  floor  tectonics  constitute  the  principal  weakness  of 
this  aspect  of  this  study.  Future  improvements  in  the  data  coverage 
should  enhance  the  reliability  of  this  analytic  tool. 

Figure  1 illustrates  the  method  of  encoding  the  grid  elements;  in 
this  case,  an  area  in  southeast  Asia.  Each  5°x5°  element  is  inspected 
for  the  visually  estimated  percentage  (to  the  nearest  10%)  of  the  en- 
closed tectonic  divisions.  In  this  example,  three  of  the  twenty  tectonic 
divisions  are  present;  shelf  sediments  40%,  Mesozoic  orogeny  30%,  and 
Post-Precambrian  undeformed  platform  30%.  Each  of  these  divisions  is 
assigned  a percentage  of  the  total  area  of  the  scene  and  the  weights  are 
entered  into  the  tectonic  matrix;  all  other  tectonic  divisions  are  as- 
signed weights  of  zero.  The  spatial  distribution  of  the  tectonic  div- 
isions within  the  grid  element  is  not  included  in  the  matrix. 

It  is  possible  that  a particular  path  through  a grid  element  will 
sample  only  a portion  of  the  enclosed  tectonics,  and  therefore  the  actual 
path  segment  velocity  may  not  be  represented  accurately  by  the  tectonic 
based  5"x5°  element  mean  velocity.  Since  it  is  likely  that  a juxtaposit- 
ion of  two  or  more  widely  dissimilar  tectonic  divisions  will  result  in  a 
mutual  reequilibration  with  time,  the  anomalous  influence  of  the  spatial 
aspect  of  the  tectonic  divisions  seems  remote,  except  when  dealing  with 
most  recent  tectono-thermal  events.  Similarly,  if  paths  which  are  much 
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greater  in  length  than  500  kilometers  are  chosen  for  analysis,  they 
will  sample  and  average  element  tectonic  influence  over  many  elements. 

The  chance  of  fortuitous  local  spatial  distributions  of  tectonic  divisions 

significantly  perturbing  results  utilizing  the  tectonic  matrix  thus 
becomes  less  likely  as  path  length  is  increased.  The  tectonic  matrix 

is  stored  on  magnetic  tape  for  easy  amendment  and  access. 

Utilization  of  this  tectonic  matrix  to  predict  Rayleigh  wave  group 
velocities  requires  the  assignment  of  a group  velocity  dispersion  char- 
acteristic to  each  of  the  tectonic  divisions  used  in  the  encoding.  One 
such  method  is  to  determine  the  characteristic  dispersion  empirically 
for  a given  area  showing  only  one  tectonic  division  and  assume  that 
the  characteristic  is  relatively  invariant. 

Although  this  research  is  an  examination  of  Rayleigh  group  vel- 
ocities across  ocean  basins,  many  of  the  seismograph  stations  are  loc- 
ated at  continental  margins.  Assignment  of  continental  group  velocities 

is  required,  therefore,  to  suppress  anomalous  edge  effects  and  unreal- 
istic velocity  gradients  at  these  boundaries.  Thus,  the  assignment  of 
'average'  continental  Rayleigh  group  velocities  to  all  twelve  continental 
divisions  is  unreasonable.  The  velocity  transition  from  ocean  to  con- 
tinent is  significantly  different  for  California  versus  New  Jersey,  for 
example.  The  group  velocities  for  the  twelve  continental  tectonic  div- 
isions are  derived  from  published  data.  All  phase  velocity  data  have  been 

converted  to  group  velocities  using  the  standard  formulation. 

Frequently,  group  velocities  assigned  to  a particular  tectonic  div- 
ision are  the  result  of  several  independent  investigations  throughout 
the  world.  Remarkably,  the  variations  in  Rayleigh  wave  group  velocities 
for  a given  tectonic  division  generally  differ  by  less  than  0.02  km/sec 
for  the  period  range  20  through  100  seconds.  Although  additional  analy- 


6 


ses  may  result  in  twelve  distinct  dispersion  characteristics  for  these 
tectonic  divisions,  only  six  could  be  differentiated  from  the  current 
literature.  The  Archaen  and  Proterozoic  shields  were  assigned  the  same 
velocities  for  the  period  range  of  20  through  100  seconds  based  on  the 
analyses  of  Brune  and  Dorman  (1963)  for  the  Canadian  Shield,  Bolt  and 
Niazi  (1964)  for  the  Australian  shields,  Bhattacharya  (1974)  for  the 
Indian  shield  and  Teng  (1974)  for  the  Asian  shields.  Similarly,  the 

Precambrian  and  post-Precambrian  platforms  and  the  shelf  sediments  were 

indistinguishable  (Brilliant,  1954;  Bolt  and  Niazi,  1964;  Ocal , 1965; 

Knopoff  and  Mai,  1967;  Biswas  and  Knopoff,  1974;  Teng,  1974;  Goncz, 

Hales  and  Muirhead,  1975;  Goncz  and  Cleary,  1976).  Mesozoic  orogenic 

zones  were  assigned  Rayleigh  wave  group  velocities  from  the  combined 
studies  of  Alexander  (1963),  James  (1971),  Knopoff  (1972),  Biswas  and 

Knopoff  (1974),  Fix  (1975),  and  North  (1976).  The  intermontane  basin 

fill  division  velocities  were  derived  from  three  studies  of  the  Basin 
and  Range  province  of  the  Western  United  States  (Biswas  and  Knopoff, 

1974;  Braile,  1975;  Keller,  1975).  Finally,  no  data  for  the  divisions 

early  Paleozoic  orogeny,  late  Paleozoic  orogeny,  and  Mesozoic  volcanics 

could  be  found  in  the  literature.  Accordingly,  these  divisions  were 

assigned  mean  continental  group  velocities  from  Oliver  (1962).  A sum- 
mary of  the  assigned  continental  group  velocities  of  fundamental  mode 
Rayleigh  waves  for  the  period  range  20  to  100  seconds  is  given  in  table 
1. 

Data  used  to  construct  and  to  test  the  functional  relationship 
between  sea  floor  age  and  Rayleigh  group  velocities  were  selected  from 
long  period  vertical  component  seismograms  of  WWSSN  and  HGLP  stations. 
All  seismograms  were  digitized  at  0.56  second  increments  with  an  ampli- 
tude resolution  of  .06  millimeters.  The  length  of  digitization  is  dep- 
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endent  on  the  dispersive  character  of  the  original  record  and  generally 
is  three  times  the  length  of  the  included  Rayleigh  train.  This  insures 
an  adequate  data  length  for  the  moving  window  analysis  (Landisman,  Dzie- 
wonski,  and  Sato,  1969)  to  be  performed  properly  and  prevents  partial 
loss  of  visual ly-obscure  low  amplitude  group  arrivals. 

Energy  levels  are  mapped  in  a group  velocity  versus  logaritnic  per- 
iod space  using  the  moving  window  analysis  of  Landisman,  Dziewonski,  and 
Sato  (1969).  The  arrival  time  of  the  peak  energy  level  of  a wave  packet 
gives  the  group  velocity  at  that  period  when  the  data  have  been  corrected 
for  instrumental  group  delay.  All  records  were  analyzed  with  a cosine 
squared  window  having  width  four  times  the  period  of  interest. 

Uncertainty  in  determining  period-specific  energy  maxima  is  depend- 
ent on  the  quality  of  the  original  seismogram  and  resolution  of  the  cosine 
squared  window.  Generally  the  velocity  uncertainty  is  approximately  0.02 
km/sec  in  the  best  cases,  but  may  be  as  large  as  0.10  km/sec  at  very 
poorly  resolved  long  periods  (greater  than  80  sec).  All  dispersion 
analyses  have  been  corrected  for  instrumental  group  delay  using  the  for- 
mulation of  Hagiwara  (1958)  with  corrections  by  Brune  (1962). 

Other  sources  of  error  in  group  velocity  interpretation  can  be  at- 
tributed to  the  seismic  source  and  transmission  path.  The  magnitudes  of 
the  source  errors  are  generally  inversely  proportional  to  the  epi- 

central  distance,  thus  by  using  teleseismic  events  at  distances  greater 
than  4000  kilometers,  most  source- related  errors  become  insignificantly 

small . 

Epicenter  mislocation  and/or  incorrect  assignment  of  origin  time 
are  potentially  the  most  significant  sources  of  error  to  a velocity  de- 
termination. Utilizing  larger  magnitude  events  (greater  than  6.0) 
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and  several  earthquakes  from  the  same  location  tend  to  reduce  this  un- 
certainty, however.  Mislocation  errors  are  estimated  to  yield  group 
velocity  uncertainties  of  the  order  of  0.02  km/sec.  This  estimate 
should  be  quite  liberal  for  most  larger  magnitude  earthquakes,  but  may 
be  conservative  for  regions  like  the  Atlantic  Indian  Rise  where  azi- 
muthal control  of  epicenter  location  is  limited. 

Group  velocity  uncertainty  attributable  to  misinterpretation  of 
initial  phase  from  the  focal  mechanism  has  been  demonstrated  by 
Knopoff  and  Schwab  (1968)  to  be  insignificant  for  large  epicentral 
distances  and  the  period  range  discussed  in  this  paper,  thus  it  is  not 
considered  as  a relevant  source  of  error.  Similarly,  other  source 
parameters  such  as  fault  length,  rupture  velocity,  and  stress  drop, 
while  affecting  the  wave-frequency  amplitudes  are  not  considered  to 
affect  significantly  the  observed  group  velocities.  Although  spceific 
source  parameters,  such  as  depth  and  rupture  velocity  can  yield  far- 
field  surface  wave  spectral  'holes'  (Turnbull,  1976).  these  energy  dis- 
continuities are  unrelated  to  path-associated  dispersion  and  can  be  dis- 
regarded. 

Discrimination  of  tectonic-related  group  velocity  variations 
through  pcean  basins  presupposes  that  other  path  related  velocity  per- 
turbations can  be  isolated  and  removed  from  the  observed  data.  The 
effects  of  horizontal  refraction  and  ocean  depth,  in  fact,  may  obscure 
completely  the  actual  group  velocity  variation  due  to  sea  floor  evolu- 
tion. A basic  assumption  necessary  for  the  computation  of  the  group 
velocities  is  that  the  surface  waves  travel  a minimum  distance  (great 

circle)  path.  This  assumption  is  strictly  correct  only  when  the  medium 
is  laterally  homogeneous  or  when  the  great  circle  path  is  normal  to 
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the  inhomogeneity  boundaries.  In  a medium  with  lateral  structural 
variations,  however,  it  is  possible  that  the  least-time  rather  than 
least-distance  path  is  responsible  for  the  arrival  on  a particular 
seismogram.  Generally  the  difference  between  the  minimum-time  and 
minimum-distance  paths  should  be  a second  order  effect  for  group  vel- 
ocities, since  the  higher  velocity  of  a refracted  path  tends  to  be 
counterbalanced  by  the  increased  path  length.  When  the  path  difference 
is  significant,  horizontal  refraction  effects  are  manifested  as 
multiple  energy  arrivals  (Capon,  1970,1971),  beating  (Pilant  and 
Knopoff,  1964),  or  unusual  amplitude  fluctuations  (McGarr,  969). 

A Significant  difference  between  a minimum-time  and  minimum- 
distance  path  occurs  when  the  velocity  gradient  normal  to  the  pro- 
pagation path  is  large.  High  crustal  and  mantle  velocity  gradients 
can  occur  in  ocean  basins  normal  to  ridge  axes  (Forsyth,  1975),  prox- 
imal to  unusual  thermal  perturbations  (hot  spots)  (McGarr,  1969)  and 
at  continental  margins  (Capon,  1970). 

The  gradient  normal  to  a ridge  axis  is  at  a maximum  during  the 
first  fifteen  million  years  of  evolution  and/or  in  basins  where  the 
spreading  rate  is  low.  The  percentage  of  the  sea  floor  where  this  can 
be  a problem  is  relatively  small  but  may  account  for  some  of  the  un- 
usual Atlantic  Ocean  test  results.  The  influence  can  be  reduced  by 
selecting  only  paths  nearly  normal  to  the  ridge  axes.  If  there  is 
significant  velocity  anisotropy  associated  with  the  spreading  directions 
(Bamford  and  Crampin,  1977;  Fuchs,  1977;  Schlue  and  Knopoff,  1977; 
Forsyth,  1975),  then  this  technique  resolves  only  the  velocities  in 


that  direction  which  may  not  be  representative  of  a more  generalized 
path.  The  simplest  solution  to  this  problem  is  to  avoid  seismograms 
which  are  complicated  by  symptomatic  beating  and  unusual  amplitude 
fluctuations. 

McGarr  (1969)  similarly  has  demonstrated  that  surface  waves  which 
have  propagation  paths  near  hot  spots  commonly  exhibit  refractive 
focusing  and  defocusing  effects.  Although  this  is  less  of  a spatial 
problem  than  are  the  ridge  axes,  it  is  a potential  problem  when  the 
recording  station  is  located  on  such  a perturbation,  for  example  KIP, 
or  if  the  propagation  path  is  aligned  with  the  perturbed  zone.  In 
general  the  effects  of  these  small  regions  can  be  considered  insig- 
nificant for  all  but  the  shortest  periods  (less  than  20  sec),  if  the 
propagation  paths  are  chosen  to  avoid  them. 

Capon  (1970)  was  able  to  show  that  surface  waves  can  undergo  con- 
siderable refraction  at  continental  margins.  The  effect  was  at  a maxi- 
mum at  LASA  (Montana  array)  when  the  waves  approached  the  margin  at 
relatively  low  angles.  The  effect  is  negligule  for  coastal  stations 
when  the  angle  of  approach  is  nearly  ninety  degrees. 

In  summary,  the  effects  of  horizontal  refraction  or  departure  from 
a minimum-distance  path  are  difficult  to  evaluate.  By  selecting  seis- 
mograms which  exhibit  few  obvious  multi-path  characteristics  and  avoiding 
paths  which  are  likely  to  create  these  effects,  the  combined  horizontal 
refraction  perturbations  to  group  velocities  can  be  considered  insig- 
nificant. 

The  influence  of  the  water  layer  on  Rayleigh  wave  dispersion  through 
ocean  basins  is  perhaps  the  most  significant  correction  which  must  be  ap- 
plied to  the  raw  data.  It  could  be  argued  that  the  increase  in  ocean 
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depth  with  the  thermal  contraction  of  the  crust  and  upper  mantle  is  a 
correlative  to  the  evolution  of  the  oceanic  lithosphere  structure 
(Sclater  and  Francheteau,  1970;  Sclater,  Anderson  and  Bell^  1971). 
However,  there  is  significant  scatter  to  the  topography-age  relation- 
ship. Failure  to  correct  for  ocean  depth  along  a defined  path  tends 
to  obscure  rather  than  elucidate  the  continuous  evolutionary  aspect 
of  the  oceanic  crust  and  upper  mantle. 

Modeling  the  absolute  effect  of  increasing  ocean  depth  in  an 
evolving  crust  and  upper  mantle  setting  would  be  the  preferred  con- 
dition. Without  a firm  understanding  of  the  exact  structural  evolu- 
tion, however,  this  model  is  not  tractable.  Assuming  a constant  crust 
upper  mantle  structure  and  varying  only  the  thickness  of  the  water 
layer  does  permit  an  evaluation  of  the  relative  influence  of  the 
water  layer  depth  on  Rayleigh  phase  and  group  velocity  dispersion. 

Since  the  evolution  of  the  sea  floor  increases  the  observed  velocities, 
while  the  increase  in  ocean  depth  decreases  the  observed  velocities, 
by  assuming  a constant  crust  upper  mantle  structure,  the  computed  in- 
fluence of  the  water  layer  will  be  an  overestimation  for  older  ocean 
basin  regions.  To  illustrate  the  effect  of  the  ocean  depth,  fundamental 
mode  Rayleigh  phase  and  group  velocities  were  calculated  for  seven  ocean 
models.  The  basic  crust  and  mantle  structure  used  throughout  the  cal- 
culations is  given  in  table  2.  The  model  differs  only  in  the  thickness 
of  the  first  layer  which  was  varied  steowise  from  zero  through  six  kilo- 
meters. The  method  of  calculating  phase  and  group  velocity  data  is 
that  of  Harkrider  (1964).  The  results  of  these  calculations  are  given 
in  table  3.  The  influence  of  the  ocean  layer  on  the  phase  and  group 
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velocities  increases  as  the  period  decreases.  Similarly,  since  the 
group  velocities  are  more  sensitive  to  structural  changes  than  are  the 
phase  velocities,  th'^y  are  significantly  more  perturbed  than  are  the 
phase  velocities  at  the  corresponding  periods.  The  group  velocity 
correction  becomes  extremely  non-linear  as  the  period  decreases,  re- 
quiring a correction  as  large  as  28%  for  deep  oceans  at  periods  of  20 
seconds.  Generally  the  correction  is  less  than  5%  or  0.02  km/ sec, 
but  this  is  significant  when  the  variation  due  to  tectonics  is  of  the 
order  of  10%.  All  observed  Rayleigh  group  velocity  dispersions  used 
in  this  study  have  been  corrected  for  mean  great-circular  path  ocean 
depth. 

In  summary,  Rayleigh  group  velocity  data,  corrected  for  instru- 
mental response  and  mean  ocean  depth,  have  a mean  RMS  residual  uncer- 
tainty of  the  order  of  0.03  km/sec.  The  most  significant  error  will 
be  introduced  by  mislocation  of  the  epicenter  and/or  incorrect  assign- 
ment of  origin  time. 

DETERMINATION  OF  THE  GROUP  VELOCITY-AGE  RELATIONSHIP 


Rayleigh  group  velocity  data  used  to  create  the  group  velocity 
versus  age  functions  were  derived  from  a combination  of  published  dis- 
persion data  and  some  new  paths.  Dispersion  data  from  the  literature 
could  be  used  only  if  they  satisfied  the  following  criteria:  (1)  the 
least-distance  path  had  to  be  chronologically  identifiable  in  terms  of 
magnetic  anomalies  and/or  deep  sea  drilling,  and  limited  to  a small 
deviation  about  the  mean  age  to  insure  'pure-path'  properties;  (2)  the 
seismograms  had  to  be  relatively  free  of  significant  apparent  multi- 
path  characteristics;  (3)  the  dispersion  curves  must  be  computed  by  a 
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moving  window-type  analysis  which  is  clearly  identifiable  as  instru- 
mentally  delay  corrected;  and  (4)  the  pure  path  data  created  by  matrix 
least-squares  reduction  of  mixed  path  data,  such  as  those  by  Forsyth 
(1975), had  to  have  demonstrably  stable  solutions  with  clearly  stated 
error  limits.  A surprisingly  small  volume  of  data  is  available  in  the 
literature  which  meet  these  criteria.  Table  4 lists  the  published  dis- 
persion analyses  used  in  this  study.  Supplemental  new  analyses  were 
done  to  spot  check  the  reliability  of  the  published  results  and  to  fill 
some  age  gaps.  All  data  compiled  were  corrected  for  instrumental  and 
ocean  depth  group  delay. 

Even  though  nearly  all  group  velocity  data  used  in  the  construct- 
ion of  the  age-velocity  functions  are  from  the  Pacific  Ocean,  where  there 
is  a reasonable  distribution  of  earthquakes  and  observatories  and  re- 
latively wide  isochronological  sea  floor  segments,  determination  of  mean 
velocities  for  a given  sea  floor  age  is  difficult  because  of  sampling 
bias  in  an  anisotropic  medium.  The  degree  of  seismic  surface  wave  vel- 
ocity anisotropy  in  the  Pacific  has  been  estimated  to  be  between  two 
and  ten  percent  by  several  authors  (Bamford  and  Crampin,  1977;  Fuchs, 
1977;  Schlue  and  Knopoff,  1977,  Forsyth,  1975).  Velocity  maxima  are 
parallel  with  spreading  directions  and  therefore  they  are  normal  to 
sea  floor  isochrons.  Conversely,  velocity  minima  are  parallel  with 
ridge  axes  and  therefore  are  parallel  to  sea  floor  isochrons. 

Unless  an  islchronological  segment  of  sea  floor  is  sampled  at  all 
azimuths,  therefore,  an  apparent  mean  velocity  will  be  biased  by  the 
particular  azimuth-dependent  velocities  used  in  the  mean  computation. 
Forsyth  (1975)  determined  mean  surface  wave  velocities  for  temporal- 
restricted  sea  floor  segments  less  than  30  million  years  old  which  are 


relatively  free  of  anisotropy  sampling  bias.  These  mean  results  from 
the  Nazca  plate  are  used  for  the  initial  points  on  the  age-velocity 
functions.  All  other  data  points  are  subject  to  aniostropy  sampling 
bias  which  cannot  be  avoided  for  the  majority  of  the  Pacific  basins. 

Assuming  that  the  surface  wave  velocity  anisotropy  of  2%  (Forsyth, 
1975)  is  representative,  and  that  the  anisotropy  is  restricted  th  the 
low  velocity  zone  (LVZ)  as  suggested  by  Schlue  and  Knopoff  (1977),  it  is 
possible  to  predict  how  specific  observations  may  be  biased.  For  ex- 
ample, the  apparent  mean  surface  wave  velocities  for  sea  floor  (Nazca) 
older  than  30  million  years  determined  by  Forsyth  (1975)  should  be 
biased  toward  higher  than  true  mean  velocities  because  sampling  must 
be  restricted  to  paths  nearly  normal  to  the  South  American  coastline 
to  avoid  significant  horizontal  refraction  effects.  Similarly,  'pure 
path'  observations  in  the  older  portions  of  the  Pacific  (greater  than 
90  million  years)  are  likely  to  be  biased  toward  lower  than  true  mean 
velocities  because  of  the  selection  of  paths  which  are  sub-parallel  to 
sea  floor  isochrons.  Without  absolute  determination  of  the  anisotropy 
correction  necessary,  however,  it  is  not  possible  to  make  this  correction 
to  the  data  in  advance. 

Construction  of  the  functional  relationships  between  sea  floor  age 
and  period-specific  apparent  mean  group  velocities  was  accomplished  by 
performing  stepwise  polynomial  regressions  on  the  data  from  table  4 at 
periods  from  30  to  100  seconds  in  10  second  increments.  The  results  of 
these  analyses  are  given  in  figure  2.  Although  the  polynomial  regres- 
sions were  able  to  fit  the  data  with  a high  degree  of  precision,  the 
physical  sense  of  these  functions  is  questionable.  The  general  trend  of 
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these  functions  is  a smooth  'logarithmic-type'  curve  in  which  velocity 
increases  with  increase  in  age  (30,50,60,70  second  curves).  The  cause 
of  the  40  second  period  function  aberrancy  is  unknown,  but  presumed  to 
be  a sampling  problem  since  neither  the  30  nor  the  50  second  period 
functions  behave  in  this  erratic  fashion.  Although  the  excursions  are 
well  within  the  bounds  of  a 1-2%  anisotropy  deviation,  preferential 
sampling  of  the  anisotropic  medium  is  not  believed  to  be  responsible 
for  the  fluctuations  in  the  40  second  curve.  On  the  other  hand,  the 
1.3%  reduction  in  velocities  of  the  80,  90,  and  100  second  period 
curves  in  the  crustal  age  region  of  90  through  120  million  years  is 
believed  to  be  caused  by  anisotropy  sampling  bias.  All  paths  sampling 
sea  floor  in  this  age  bracket  are  subparallel  to  the  sea  floor  iso- 
chrons. 

The  temporal  limits  of  the  seven  oceanic  tectonic  divisions  (other 
than  island  arcs)  which  are  used  in  this  study  are  indicated  by  the 
vertical  lines  on  the  age  versus  apparent  mean  group  velocity  functions 
of  figure  2.  A set  of  mean  apparent  group  velocities  for  periods  from 
30  through  100  seconds  is  assigned  to  each  tectonic  division  by  averag- 
ing the  apparent  velocities  within  each  division  limit.  The  assigned 
apparent  mean  group  velocities  for  each  division  are  given  in  table  1. 
Although  there  is  a set  of  velocities  for  a 20  second  period,  this  func- 
tion was  too  erratic  to  be  considered  meaningful. 

THE  GLOBAL  RAYLEIGH  GROUP  VELOCITY  MATRIX 

Apparent  mean  Rayleigh  group  velocities  through  each  5°x5°  geo- 
graphic block  are  derived  by  summing  the  products  of  the  tectonic  ele- 
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ment  percentages  and  the  period-specific  mean  group  velocities  assigned 
to  each  tectonic  element;  or  in  matrix  notation: 

where  V is  the  period-specific  mean  group  velocity  matrix  with  dimen- 
sions 8x36x72,  T is  the  tectonic  matrix  with  dimensions  20x36x72  and  V 
is  the  period-specific,  tectonic-specific  mean  group  velocity  matrix 
with  dimensions  8x20. 

Eight  Mercator  conformal  projections  of  the  contoured  mean  Rayleigh 
wave  group  velocities  for  the  periods  30  through  100  seconds  were  pro- 
duced from  the  velocity  matrix.  Figures  3,4,  and  5 illustrate  the  velo- 
city contouring  for  40,  70,  and  100  second  Rayleigh  waves  respectively 
(13"x26"  blackline  reproductions  of  all  eight  maps  are  available  from 
the  author  on  request).  These  maps  are  a convenient  data  display  pro- 
viding both  a base  for  nonautomatic  'ray  tracing'  testing  of  the  model 
against  observed  data.  All  maps  are  contoured  at  an  interval  of  0.05 
km/sec.  Each  map  displays  the  velocity  data  from  75°  north  latitude 
to  75°  south  latitude.  There  is  a duplication  of  65°  east  longitude 
through  105°  east  longitude  on  the  left  and  right  hand  margins  of  the 
maps.  East  longitudes  0°,  180°,  and  270°  and  colatitudes  45°,  90°  and 
135°  provide  a reference  framework  on  all  maps. 

TESTS  OF  THE  RAYLEIGH  GROUP  VELOCITY  PREDICTIVE  MODEL 


PACIFIC  OCEAN  BASINS 
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Knopoff  (1969)  demonstrated  that  group  travel  times  in  inhomo- 
geneous dispersive  media  are  equivalent  to  the  sum  of  the  group  vel- 
ocity delays  within  each  homogeneous  path  segment.  For  this  study, 
path  segments  limited  by  mapped  successive  velocity  contours  are 
defined  as  homogeneous  divisions  and  assigned  velocities  equivalent 
to  the  mean  of  the  bounding  contours.  Paths  which  cross  no  mapped 
velocity  contours  are  assigned  velocities  dependent  on  the  spatial 
configuration  of  the  path  within  the  contoured  zone.  Contoured  areas 
with  widths  less  than  1500  kilometers  are  assigned  velocities  equi- 
valent to  the  mean  of  the  bounding  contours.  On  the  other  hand, 
contoured  areas  with  widths  exceeding  1500  kilometers,  such  as  those 
of  the  western  Pacific,  are  roughly  divided  into  three  subregions. 

The  two  marginal  segments  are  assigned  velocities  equivalent  to  the 
nearest  contour;  the  central  segment  is  assigned  a mean  velocity.  Thus 
for  example,  a path  subparallel  to  a 4.0  km/sec  contour  located  in 
a wide  ’iso-velocity'  region,  and  having  the  majority  of  its  length 
within  the  pne  third  proximal  region  would  be  assigned  a velocity  of 
4.0  km/sec  rather  than  a mean  velocity. 

In  addition,  implicit  in  these  analyses  is  the  assumption  that 
all  paths  are  least-distance,  and  that  the  composite  group  velocity 
along  a given  path  is  represented  by  the  sum  of  the  assigned  segment 
group  delays.  Thus,  by  overlaying  each  of  the  eight  Rayleigh  velocity 
contour  maps  with  correctly  scaled  Mercator  conformal  projections  of 
least-distance  paths,  discrete  period,  path  related  dispersion  curves 
are  computed  by  summing  the  product  of  path  segment  group  velocities 
and  corresponding  path  length  percentages  over  the  total  path  length. 
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These  predicted  dispersion  curves,  corrected  for  mean  bathymetry,  are 
directly  comparable  with  the  instrumental -delay-corrected  observed 
group  velocities  which  are  determined  for  the  same  path  by  moving  window 
analysis. 

Seventeen  random  paths  were  selected  for  the  Pacific  Ocean  basins. 
These  paths,  illustrated  in  figure  6 , collectively  sample  the  Pacific 
basins  in  a non-preferential  manner  with  respect  to  tectonics,  possible 
anisotropy,  and  bathymetry.  Rayleigh  group  velocities  are  predicted  for 
each  of  these  paths  by  the  previously  discussed  technique.  These  pre- 
dicted dispersion  curves  are  compared  with  moving  window  analyses  of 
path-equivalent  observed  dispersion  curves  in  figures  7 through  10. 

The  six  dispersion  analyses  of  figure  7 are  presented  in  a format 
unlike  the  subsequent  comparative  diagrams,  to  illustrate  not  only  the 
comparison  of  this  model  and  observed  data,  but  also  the  comparison  of 
this  tectonic-based  model  with  three  other  standard  ocean  models.  Ob- 
served data,  corrected  for  bathymetry  effects,  are  illustrated  by  the 
closed  circles;  uncorrected  for  bathymetry  effects  by  the  open  circles. 
Solid  lines  represent  dispersion  analyses  predicted  by  this  model;  dot 
and  dash  lines  represent  dispersion  using  the  Harkrider-Anderson  ocean 
model  (1964);  dashed  lines  represent  dispersion  using  standard  ocean 
model  CIT  llA  (Kanamori  and  Abe,  1968);  and  dotted  lines  represent  dis- 
persion using  standard  ocean  model  8099  (Kanamori  and  Abe,  1968).  Only 
in  figure  7 is  the  bathymetry  correction  applied  to  the  observed  data. 
All  other  diagrams  have  the  bathymetry  correction  applied  to  the  pre- 
dicted data  to  correctly  locate  it  in  the  observed  period  versus  group 
velocity  space. 
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Figure  7 clearly  demonstrates  the  inadequacy  of  previous  fixed- 
structure  models  to  fit  the  Rayleigh  group  velocity  data.  Although 
the  Harkrider-Anderson  model  may  fit  Rayleigh  group  velocity  data  for 
old  ocean  basins  (Fig. 7, KIP  1)  with  relatively  small  residuals,  the 
model  fails  significantly  for  a younger  oceanic  region  (Fig.  7,  PEL  6). 
Similarly,  fixed  structure  models  like  CIT  llA  and  8099  may  yield  small 
residuals  over  a very  limited  spectral  band,  but  fail  to  fit  the  ob- 
served data  over  a broader  spectrum.  On  the  other  hand,  the  tectonic- 
based  model  has  an  implicitly  evolving  structure  with  time,  albeit  not 
explicitly  stated,  which  permits  it  to  fit  both  old  and  young  ocean 
crust  and  mantle  Rayleigh  group  velocities  over  a wide  spectrum.  This 
is  clearly  apparent  in  figure  7. 

The  predictive  dispersion  model  significantly  overestimates  the 
group  velocities  for  the  path  PEL  6 (Fig.  7).  Interestingly,  however,  a 
significant  portion  of  this  path  traverses  the  Nazca  ridge  which  has 
been  shown  recently  by  Cutler  (1977)  to  have  an  anomalously  thick  crust 
with  modeled  low  density  core  extending  to  at  least  20  kilometers. 

Thus,  the  apparent  failure  of  the  tectonic-based  model  to  predict  ac- 
curately the  dispersion  along  this  path  independently  illustrates  the 
Nazca  Ridge  to  be  an  anomalous  structure  even  into  the  uppermost  mantle. 
Because  the  predictive  model  is  dependent  on  surficial  tectonics,  sig- 
nificant overestimations  of  observed  dispersed  Rayleigh  group  velocities 
may  suggest  locations  of  paleo  mantle  upwellings  which  no  longer  have 
geophysical  expression  at  the  surface.  This  hypothesis  is  supported 
more  convincingly  if  the  anomaly  extends  to  longer  periods  (greater  than 
40  sec)  suggesting  a deeper  seated  cause  of  the  negative  velocity  res- 
iduals. 
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I The  remaining  dispersion  analyses  in  this  paper  compare  the  observed 

[ data  to  the  tectonic  based  model  only.  The  observed  data  are  presented 

I as  energy  contours  in  an  instrumental  group  delay  corrected  group  velo- 

city versus  logarithmic-period  space.  Although  the  entire  region  is  cap- 
able of  being  contoured,  only  the  most  significant  energy  regions  are 

t contoured  to  reduce  unnecessary  cluttering  of  the  diagrams.  Depth-cor- 

1- 

f rected  predicted-mean  group  velocity  data  for  the  same  great-circular 

path  are  presented  as  solid  circles.  The  diameter  of  the  circle  is 
approximately  half  the  RMS  uncertainty  of  0.03  km/sec  for  the  determin- 
ations. The  vertical  bars  extending  from  the  mean  velocities  correspond 
to  a 2%  velocity  uncertainty  at  4.0  km/sec  (0.08  km/sec).  If  the 
degree  of  surface  wave  velocity  anisotropy  in  ocean  basins  does  not  ex- 
ceed 2%,  as  suggested  by  Forsyth  (1975),  then  these  error  bars  should 
encompass  the  possible  anisotropy-related  scatter.  Figures  8 through  10 
illustrate  the  remaining  Pacific  Ocean  tests. 

Generally  the  correspondence  of  the  predicted  and  observed  dis- 
persion data  is  remarkably  high.  A systematic  overestimation  of  the 
30  and  frequently  40  second  period  group  velocities  suggests  that  the 
depth  correction  by  mean  ocean  depth  may  not  be  adequate.  Alternatively 
this  could  be  a problem  with  the  model,  but  without  performing  a step- 
wise depth  correction,  the  source  of  this  error  is  unknown.  Dispersion 
curves  for  paths  crossing  the  Philippine  plate  are  illustrated  in  figure 
10-B,  10-C,  and  10-D.  Two  of  these  paths  cross  the  width  of  the  plate, 
the  third  crosses  the  length  (See  figure  6.).  Although  a predicted 
dispersion  is  derived  for  these  paths,  all  three  moving  window  analyses 
of  the  observed  data  are  severely  distorted.  Energy  contours  at  periods 
less  than  40  seconds  indicate  the  presence,  of  'spectral  holes'  and  the 


resolving  kernals  at  periods  greater  than  40  seconds  are  extremely 
broad,  indicating  a severe  attenuation  of  longer  period  waves.  Studies 
by  deBremaecker  (1958),  Hudson  and  Knopoff  (1964),  Pilant,  Knopoff,  and 
Schwab  (1964),  and  Mai  and  Knopoff  (1966)  have  shown  that  transmission 
and  reflection  of  Rayleigh  waves  at  abrupt  velocity  contrasts,  especial- 
ly at  corners  can  cause  severe  distortion  of  observed  Rayleigh  waves. 
Considering  that  MAT  and  GUA  are  both  located  at  critical  points  for 
effective  Rayleigh  wave  transmission  (deBremaecker,  1958),  the  lack  of 
resolution  on  the  Philippine  plate  tests  is  understandable. 

INDIAN  OCEAN  BASIN 


Four  paths  crossing  the  Indian  Ocean  are  illustrated  in  figure  11. 
The  associated  dispersion  analyses  are  illustrated  in  figure  12.  Only 
the  Rayleigh  group  velocity  dispersion  recorded  at  MUN  (figure  12-D) 
has  consistently  small  velocity  residuals  with  respect  to  the  model  pre- 
diction. All  other  analyses  consistently  predict  group  velocities  ap- 
proximately 0.10  km/sec  less  than  those  actually  observed.  Although 
a systematic  regional  mislocation  error  cannot  be  entirely  ruled  out 
as  the  cause  of  this  anomaly,  the  simultaneous  recording  of  the  same 
earthquake  at  CHG  (fig.  12-C)  and  MUN  (fig.  12-D)  suggests  that  the 
error  is  path  related  and  not  the  result  of  mislocation.  Since  the  paths 
to  CHG,  KOD,  and  LEM  (figure  12.)  sample  relatively  the  same  aged 
material  (Johnson,  Powell  and  Veevers,  1976),  it  is  possible  that  the 
relatively  consistent  positive  velocity  offset  is  a result  of  incorrect 
age  assignment.  This  is  a relatively  old  section  of  sea  floor,  how- 
ever, and  the  misassignment  of  age  would  have  to  be  extreme  to  account 
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for  the  magnitude  if  the  error.  Similarly,  the  magnitude  of  this  error 
is  too  large  to  be  the  result  of  an  incorrect  depth  correction.  Further 
discussion  of  this  positive  anomaly  is  deferred  to  the  discussion  and 
speculation  section. 


ATLANTIC  OCEAN  BASINS 

Six  Atlantic  Ocean  test  paths  are  illustrated  in  the  Mercator 
conformal  projection  of  figure  13.  The  dispersion  analyses  for  the 
three  South  Atlantic  basins  tests  are  illustrated  in  figures  14A,  14B, 
and  14C.  The  dispersion  analysis  ,14A,  is  for  a path  traversing  the 
Cape  Basin  of  the  South  Atlantic  Ocean,  and  is  characterized  by  pos- 
itive observed-minus-predicted  velocity  residuals  approximately  equi- 
valent to  those  observed  for  the  Indian  Ocean.  Both  dispersion  curves 
14B  and  14C  are  for  paths  either  partly  (14B)  or  entirely  (14C)  con- 
fined to  the  Brazil  Basin.  Both  dispersion  curves  exhibit  extremely 
large  positive  residuals  with  respect  to  the  model;  14C  larger  than 
14B.  Since  P-wave  travel  times  for  both  events  have  essentially  no 
residuals,  timing  and/or  mislocation  errors  can  be  discounted.  Sim- 
ilarly, there  are  no  input  errors  in  the  moving  window  analyses,  thus 
these  anomalies  appear  to  be  real  and  related  to  the  Brazil  Basin. 

North  Atlantic  tests  are  illustrated  in  figures  14D,  15A,  and 
15B.  Inadequate  resolution  of  the  moving  window  analysis  for  the 
recording  of  the  earthquake  at  TRN  (fig.  15A)  appear  to  be  character- 
istic for  that  station  since  similar  observations  have  been  made  by 
other  authors  (Weidner,  1972).  Since  the  resolution  difficulties  pre- 
clude interpretation  of  TRN  recordings  of  mid-Atlantic  earthquakes. 
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this  station  is  unsuitable  for  surface  wave  analyses  of  Atlantic  paths. 
The  other  two  dispersion  curves  for  the  North  Atlantic,  14D  and  15B, 
are  not  confined  to  specific  basins  but  have  paths  subnormal  (14D)  and 
subparallel  (15B)  to  the  mid-Atlantic  ridge.  Both  of  these  paths  have 
characteristic  velocity  residuals  (0.1-0. 2 km/sec)  larger  than  those 
observed  for  the  Cape  Basin  (0.05-0.1  km/ sec)  in  the  South  Atlantic  or 
the  Indian  Ocean  (0.05-0.1  km/sec),  but  significantly  less  than  the 


residuals  observed  for  the  Brazil  Basin. 

In  summary,  tests  of  the  tectonic-based  Rayleigh  wave  group  vel- 
ocity predictor  have  been  conducted  for  many  of  the  major  ocean  basins. 
Only  in  the  Pacific  basins  does  the  model  consistently  predict  the 
dispersion  curves  for  given  paths  with  extremely  small  errors.  In 
all  other  ocean  basins,  the  model  systematically  underestimates  the 
observed  velocities. 

DISCUSSIONS  AND  SPECULATIONS 

The  relative  success  of  this  Rayleigh  group  velocity  model  clearly 
indicates  that  the  first  order  variations  in  observed  oceanic  disper- 
sion are  related  to  sea  floor  age.  In  the  Pacific  Ocean,  there  is  a 
10%  increase  in  the  group  velocity  of  Rayleigh  waves  with  a period  of 
40  seconds  over  the  span  of  140  million  years.  Approximately  6%  of 
this  10%  increase  occurs  within  the  first  30  million  years  of  evolution. 
In  comparison,  there  is  a 7%  increase  in  group  velocity  of  Rayleigh 
waves  with  a period  of  90  seconds  over  the  span  of  140  million  years, 
with  only  3%  of  the  7%  increase  occurring  within  the  first  30  million 
years. 
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It  is  gratifying  that  the  Rayleigh  predictor  model  is  as 
successful  as  it  appears  to  be  in  the  Pacific  basins.  However,  the 
failures  of  the  model  pose  intriguing  problems. 

These  can  be 

separated  into  explanations  of  negative  residuals  and  explanations 
of  positive  residuals. 

The  one  path  displaying  systematic  negative  observed-minus-pre- 
dicted Rayleigh  group  velocity  residuals  for  the  Pacific  Ocean  tra- 
verses the  Nazca  Ridge  which  has  been  shown  to  have  an  anomalous 
crustal  structure  and  low  density  core  (Cutler,  1977).  This  suggests 
a previously  unexpected  application  of  the  Rayleigh  predictor  model, 
i.e.  the  search  for  'fossil'  spreading  regions.  These  ancient  thermal 
perturbations  should  retain  a mantle  signature  significantly  beyond 
the  time  that  surface  expression  (topography)  of  the  divergent  zone 
has  disappeared  (approximately  35  million  years).  Thus,  for  example, 
it  may  be  possible  to  located  the  mantle  'root'  of  the  Kula  Ridge 
even  though  the  surface  geophysical  signature  (magnetic  anomaly  pat- 
terns) of  the  axis  of  this  spreading  center  has  been  spatially  decoupled. 
The  location  of  fossil  upwelling  zones  in  the  mantle,  which  should  be 
fixed  spatially  for  a longer  period  of  time  than  their  surface  cor- 
relative features,  can  offer  additional  constraints  to  possible  paleo- 
plate  reconstructions.  This  is  especially  true  when  episodes  of  plate 
motion  are  characterized  by  distinct  directions  of  motion  which  are 
time  dependent,  and  subsequently  deactivated  ridges  are  spatially 
migrated  as  has  occurred  in  the  Pacific. 

The  systematic  positive  velocity  residuals  with  respect  to  the 
Rayleigh  predictor  which  occur  in  all  ocean  basins  except  the  Pacific 


constitute  a second  intriguing  problem.  Although  the  observed  higher 
Rayleigh  velocities  of  the  non-Pacific  basins  may  be  a function  of 
compositional  variation,  the  relative  world-wide  homogeneity  of  sam- 
pled oceanic  lithosphere  (Chung,  1976)  does  not  support  this  proposi- 
tion. Spatial  variability  of  mantle  heat  flux  (Pollack  and  Chapman, 
1977a, b)  appears  to  be  a more  likely  cause  of  the  Rayleigh  group 
velocity  variation. 

Two  alternative  explanations  for  the  lower  Atlantic  and  Indian 
Ocean  crust  and  upper  mantle  thermal  structure,  and  thus  higher  Ray- 
leigh group  velocities,  are  possible.  Either  there  may  be  a lower 
initial  thermal  input,  or  there  is  a mechanism  for  extraction  of  the 
heat  which  is  not  present  in  the  Pacific.  If  the  former  explanation 
is  correct  as  suggested  by  Pollack  and  Chapman  (1977b),  then  the  re- 
duction of  partial  melting  due  to  reduced  oceanic  geothermal  gradients 
is  a relatively  time-invariant  influence  and  provides  a constant  off- 
set of  the  group  velocity  versus  age  functions  of  figure  2.  Thus,  for 
example,  paths  through  equivalently  aged  lithosphere  in  the  Atlantic 
versus  the  Pacific  would  display  consistently  higher  velocities  in 
the  former  basins  than  the  latter,  and  regardless  of  age,  the  offset 
would  be  the  same. 

Alternatively,  the  Atlantic  and  Indian  Oceans  differ  from  the 
Pacific  in  another  important  aspect,  i.e.  the  former  basins  have  been 
locations  occupied,  perhaps  repeatedly,  by  continental  lithosphere, 
whereas  the  latter  has  apparently  remained  an  oceanic  lithosphere  realm 
throughout  its  existence.  The  cooler  cratonic  lithosphere,  proposed 
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by  Pollack  and  Chapman  (1977b)  to  be  in  excess  of  300  kilometers  thick, 
has  a profound  influence  on  the  thermal  structure  of  the  lower  crust 
and  upper  mantle.  The  degree  of  influence  should  be  proportional  to 
the  distance  from  the  spreading  center  to  the  cratonic  margin.  If  the 
proximity  to  cratonic  elements  rather  than  lower  spreading  center  temp- 
eratures is  the  cause  of  the  higher  Rayleigh  group  velocities  in  ocean  basins 
occupied  by  continents  prior  to  drifting, then  the  group  velocity  versus  age 
functions  (figure  2.)  would  be  influenced  in  a temporally  'lessening 
manner.  The  oldest  ocean  basins,  closest  to  the  initial  rifting  would  be 
most  pronouncedly  affected,  while  the  present  mid-Atlantic  ridge 
should  show  little  or  no  influence  of  era ton  proximity  on  its  geotherms. 

This  could  explain  why  the  Brazil,  and  possibly  Angola  Basins  would 
have  very  fast  Rayleigh  group  velocities  since  they  were  closest  to 
cratonic  elements  contemporaneous  with  the  breakup  of  Gondwana.  Sim- 
ilarly, the  basins  of  the  North  Atlantic  would  have  fast  Rayleigh  vel- 
ocities due  to  the  proximity  with  North  American  and  African  cratonic 
regions.  ’ However,  these  basins  should  be  slower  than  the  Brazil  Basin 
since  there  was  a younger  tectonically  disturbed  lithospheric  buffer 
between  these  cratons  and  the  incipient  spreading  center. 

Additional  testing  of  the  Atlantic  and  Indian  Ocean  basins  to 
attempt  to  differentiate  between  these  alternative  explanations  for 
positive  residuals  from  the  Rayleigh  predictor  model  is  a currently 
active  research  effort  at  the  University  of  Michigan.  In  all  likelihood,  the 
higher  velocities  characteristic  of  the  non-Pacific  ocean  basins  might 
be  the  result  of  a combination  of  lower  initial  thermal  flux  and 

thermal  diffraction  due  to  continental  proximity  influence. 
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APPENDIX  Sources  of  Tectonic  Data 

Continents 

North  America:  Tectonic  Map  of  North  America  (1969).  United 
States  Geological  Survey,  scale  1:5,000,000  compiled  by  P,  B, 
King. 

South  America;  Geologic  Map  of  South  America  (1964),  Com- 
mission of  the  Geologic  Map  of  the  World,  scale  1:5,  000,000, 
General  coordinator  A.R.  Lamego;  Tectonic  Map  of  Brazil  (1971), 
Ministry  of  Mines  and  Energy,  National  Department  of  Mineral 
Production,  Brazil,  scale  1:5,000,000,  coordinator  E.O.  Ferreira. 

Antarctica  : Geologic  Map  of  Antarctica  (1971),  American  Geo- 
graphical Society,  scale  1,  5,000,000,  compiled  by  C,  Craddock. 

Africa:  International  Tectonic  Map  of  Africa  (1968),  Association 
of  African  Geologic  Surveys  and  United  Nations  Educational  Sci- 
entific and  Cultural  Organization,  scale  1:5,000,000,  coordinator 
G.  Choubert. 

Australia  ; Tectonic  Map  of  Australia  and  New  Guinea  (1971), 
Geological  Society  of  Australia,  scale  1:5,000,000,  Sidney. 

Europe  and  As^  ; Tectonic  Map  of  Eurasia  (1966),  Geologic 
Institute  of  the  Academy  of  Sciences  of  the  U.S.S.R.,  scale  1: 
5,000,000,  chief  editor  A.L.  Yanshin. 

Oceans 

General;  Magnetic  Lineations  of  the  Oceans  (1974),  Geologic. 
Society  of  America,  compiled  by  W.  C.  Pitman  III,  R.  L.  Lar- 
son, and  E.  M.  Herron. 

Arctic  ; Herron,  E,  M.,  J.  F.  Dewey,  and  W.  C.  Pitman  III 
( 1974), 'Plate  Tectonics  Model  for  the  Evolution  of  the  Arctic', 

Geol . (2),  p.  377-380. 
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Sources  of  Tectonic  Data  continued 


Indian  : Johnson,  B.  D.,  C.  Mca.  Powell,  and  J.  J.  Veevers 
( 1976), 'Spreading  History  of  the  Eastern  Indian  Ocean  and  Greater 
India's  Northward  Flight  from  Antarctica  and  Australia',  Geol. 

Soc . Amcr . Bull.  (87),  p.  1560-1566. 

Marginal  Basins  ; Sclater,  J.G.,  U.G.  Ritter,  and  F.  S.  Dixon 
(1972), 'Heat  Flow  in  the  Southwestern  Pacific',  Jour . Geophys . Res . 
(77),  p.  5697-5704.;  Sclater,  J .G . (1972), 'Heat  flow  and  Ele- 
vation of  the  Marginal  Basins  of  the  Western  Pacific',  Jour.  Geo- 
phys . Res . (77),  p,  5705-5719.  ; Cooper,  A.K.,  M.S.  Marlow, 
and  D.W.  Scholl  (1976),'  Mesozoic  Magnetic  Lineations  in  the 
Bearing  Sea  Marginal  Basin',  Jcur.  Geophys . Res . (81),  p,  1916- 
1934.;  Lawver,  L.  A.,  J.  R.  Curray,  and  D.G.  Moore  (1975),' 
Magnetics  in  the  Adaman  Sea  and  the  Effect  of  High  Sedimentation 
Rates',  ^.O.S.  (56),  p.l064.;  Packham,  G.H.,  and  D.  A.  Falvey 
(1971), 'An  Hypothesis  for  the  Formation  of  Marginal  Seas  in  the 
Western  Pacific',  Tectonophys.  (11),  p.  79-110.;  VanderVoo,  R., 
F.J.  Mauk,  and  R.B.  French  ( 1976), 'Permian  Triassic  Contin- 
ental Configurations  and  the  Origin  of  the  Gulf  of  Mexico',  Geol  . 
(3),  p.  177-179.;  Initial  Reports  of  the  Deep  .Sea  Drilling  Project 
Volumes  Prepared  for  the  National  Science  Foundation  under 
Contract  C-482 . 
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CAPTIONS  FOR  FIGURES 

Figure  1.  5°x  5°  grid  element  tectonic  encoding  method 

Figure  2.  Apparent  Mean  Rayleigh  Group  Velocity  Versus  Sea  Floor  Age. 

Numerical  designations  on  each  curve  indicate  the  period  of 
the  wave  in  seconds.  Vertical  bars  separate  tectonic  divisions 

Figure  3.  40  Second  Period  Rayleigh  Wave  Group  Velocities.  Contour  in- 

terval is  0.05  km/sec. 

Figure  4.  70  Second  Period  Rayleigh  Wave  Group  Velocities.  Contour  in- 

terval is  0.05  km/ sec. 

Figure  5.  100  Second  Period  Rayleigh  Wave  Group  Velocities.  Contour 

interval  is  0.05  km/sec. 

Figure  6.  Mercator  Conformal  Projection  of  Pacific  Ocean  Test  Paths. 

Figure  7.  Six  Pacific  Dispersion  Analyses.  Closed  circles  indicate  data 
corrected  for  bathymetry,  open  circles  indicate  data  uncorrect- 
ed. Solid  lines  are  dispersions  predicted  from  this  model. 

Dot  and  Dash  Lines  are  dispersions  with  Harkrider-Anderson 
Ocean  Model  (1964).  Dotted  lines  are  dispersions  with  stan- 
dard ocean  model  8099,  Dashed  lines  are  dispersions  with 
standard  ocean  model  CIT  11 A. 

Figure  8.  Pacific  Ocean  Tests.  A.  KIP-12.5  S,  77.8  W,  epicentral  distance 
9514  km.;  B.  LON-8.4  S,  154.0  E,  epicentral  distance  10230  km.; 
C.  KIP-53.0  N,  171.1  E,  epicentral  distance  4384  km.;  D.  AFI- 
39.8  N,  143.1  E,  epicentral  distance  7553  km. 

Figure  9.  Pacific  Ocean  Tests.  A.  KIP-14.7  S,  175.6  E,  epicentral  dis- 
tance 4429  Km.;  B.  KIP-15.7  S,  172.64  W,  epicentral  distance 
4397  km. 
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CAPTIONS  FOR 

Figure  10. 

Figure  11. 
Figure  12. 

Figure  13. 
Figure  14. 

Figure  15. 


FIGURES  CONTINUED 


Pacific  Ocean  and  Philippine  Plate  Tests.  A.  KIP-14.8  S, 

174.8  W,  epicentral  distance  4400  km.;  B MAT-1.2  N,  126.3 
E,  epicentral  distance  4094  km.;  C.  GUA-5.1  N,  126.8  E,  epi- 
central distance  2198  km.;  D.  6UA-1.2  N,  126.3  E,  epicentral 
distance  2455  km. 

Mercator  Conformal  Projection  of  Indian  Ocean  Test  Paths. 

Indian  Ocean  Tests.  A.  KOD-47.9  S,  99.7  E,  epicentral  dis- 
tance 6796  km.;  B.  LEM-26.1  S,  84.3  E,  epicentral  distance 
3261  km.;  C.  CHG-41.7  S,  79.9  E,  epicentral  distance  6954  km. 
D.  MUN-41.7  S,  79.9  E,  epicentral  distance  3384  km. 

Mercator  Conformal  Projection  of  Atlantic  Ocean  Test  Paths. 

South  and  North  Atlantic  Tests.  A.  SDB-58.5  S,  13.6  W,  epi- 
central distance  5307  km.;  B.  NAT-54.4  S,  2.5  E,  epicentral 
distance  6416  km.;  C.  NAT-11.7  S,  14.0  W,  epicentral  distance 
2421  km.;  D.  TLO-10.6  N,  63.4  W,  epicentral  distance  6653  km. 

North  Atlantic  Tests.  A.  TRN-15.3  N,  45.9  W,  epicentral  dis- 
tance 1756  km.;  B.  AKU-00.9  S,  22.1  W,  epicentral  distance 
7384  km. 


Table  ?..  Standard  Oceanic  Crust  and  Upper  Mantle  Model 


Layer  Thickness 

V 

V 

Rho 

(kilometers) 

P 

km/sec 

s 

km/sec 

g/cm^ 

5.0 

1,52 

0.0 

1.03 

1.0 

2.10 

l.OO 

2.10 

2-5 

6.41 

3.70 

3.06 

2.5 

6.411 

3.702 

3.07 

3.0 

8.11 

4.61 

3.40 

5.5 

8.  12 

4.61 

3.40 

17.5 

8.  12 

4.61 

3.40 

20.0 

8,01 

4,56 

3.37 

20.0 

7.95 

4,56 

3.37 

20.0 

7.71 

4.40 

3.37 

20.0 

7.68 

4.34 

3.33 

20.0 

7 . 7 / 

4.34 

3.33 

20.0 

7.85 

4.34 

3.33 

20.0 

8.10 

4.45 

3.33 

29.0 

8.  12 

4.45 

3.33 

20.0 

8.  12 

4.45 

3,35 

20.0 

8.  12 

4.45 

3.36 

20.0 

8.12 

4.45 

3,37 

20.0 

8.12. 

4.45 

3.38 

20.0 

8.24 

4 . -1  5 

3.39 

10.0 

8.50 

4.53 

3.44 

20.0 

8.  56 

4,56 

3.50 

25.0 

8.75 

4.795 

3.68 

20.0 

9.  15 

5.04 

3.88 

10.0 

9.43 

5.217 

3.90 

20.0 

9.76 

5.40 

3.92 

25.0 

9.77 

5.40 

3.93 

25.0 

9.78 

5,40 

3.95 

25.0 

9.78 

5,40 

3.96 

25.0 

9.78 

5,40 

3.99 

■H 


Table  3.  Phase  and  Group  Velocities  for  the  standard  ocean  model  of  table  2 varying  the 
thickness  of  layer  1 (water  layer)  from  0 through  6 kilometers. 
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University,  New  York. 
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(1)  Pa  and  Sa  Phases  Revis ited  by  Frederick  J.  Mauk  sub- 
mitted for  the  47th  Annual  Meeting  of  the  Eastern  Section  of  the 
Scismological  Society  of  America 

"Examination  of  four  months  of  seismograms  from  the 
White  Pine  HGLP  observatory  has  revealed  tlie  occurrence 
of  vertical  mode  Pa  and  Sa  phases  for  an  azimuthally 
restricted  population  of  teleseismic  events.  These  earth- 
quakes range  in  distance  from  50  to  106  degrees  and  in 
depth  from  30  to  60  kilometers.  The  predominant  periods 
are  15  and  25-30  seconds  for  Pa  and  Sa , respectively.  The 
apparent  velocities  of  both  Pa  and  Sa  show  strong  dependence 
on  the  type  of  lithosphere  through  which  the  waves  pass.  Pa 
ranges  from  8.05  km/sec  to  8.28  km/sec  and  Sa  from  4.58  km/sec 
to  4.54  km/sec.  The  predominant  north-south  azimuths  to  the 
events  yielding  Pa  and  Sa  display  a remarkable  correlation 
with  regions  of  thick  lithosphere  proposed  by  Chapman  and 
Pollack  (1975).  This  suggests  that  the  "whisper  gallery" 
effect  of  Press  and  Ewing  (1955)  rather  than  the  low  velocity 
zone  wave  guide  of  Caloi  (1953)  is  the  propagation  mechanism." 


CHARACTERISTICS  OF  HGLP  BODY  WAVE 
SPECTRA 

B.  Radovich  Williums  (Dept,  of  Geology  and 

Mineralogy,  Univ.  of  Michigan,  Ann  Arbor, 
Michigan  48109) 

Body  wave  FFT  and  Maximum  Entropy 
spectra  recorded  on  the  11  station  HGLP  net- 
work for  several  large  events  (mi)>  6.0)  in 
1974  have  been  determined.  This  is  part  of 
a study  to  relate  the  regionally  controlled 
Mg-mb  value  to  corner  frequency,  and  to 
rise  time  estimated  from  seismic  moment 
(Kanamori,  1972),  The  FFT  spectra  for  P 
and  SH,  corrected  for  instrument  response, 
are  consistent  with  theories  that  predict  a low 
frequency  rise  and  flattening  of  spectra.  The 
PP  phase  spectra  are  very  similar  to  P imply- 
ing little  path  distortions  at  these  long  periods, 
and  suggests  PP  can  be  used  for  source  para- 
meter studies  at  A > 100".  The  position  of  the 
corner  frequency  with  the  HGLP,  FFT  spec- 
tra is  no  more  accurate  than  with  the  WWSSN 
estimates  due  to  poor  resolution  at  long  periods 
inherent  in  the  FFT  method.  The  Maximum 
Entropy  spectra  arc  also  calcvilated  for  this 
data  set  and  show  a smoother  low  frequency 
rise  than  the  FFT  data,  but  do  not  eniiam.c  the 
flattening.  This  suggests  that  long  period 
noise  may  be  covering  the  position  of  the  corner 
frequency  and  pass  band  filters  .are  being 
designed  to  suppress  this  effect. 
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SEISMIC  MOMENTS  FROM  BRUNE'S  AR  METHOD 
AND  DETERMINATIONS  FOR  INTER- INTRAPLATE 

EVENTS 

1.  032499WILLIAiMS 

Barbara  R.  VJilliams  (Department  of  Geology 

and  Mineralogy,  University  of  Michigan,  2.  Eastern  Section 
Ann  Arbor,  MI  48109)  of  SSA,  1976 


Seismic  moments  of  20  shallow  earth- 
quakes in  1974-1975  with  5.0<M3<8.0 
have  been  calculated  using  a modified 
Brune's  AR  parameter  method.  Several 
major  events  in  this  data  set  are  of  an 
Intraplate  nature,  i.e.  not  on  well  de- 
fined plate  boundaries,  including  the 
Leeward  Is.  1974  event  with  new  Mq  =2,2 
xl026  dyne-cm,  the  Utah-Idaho  1975 
event  with  Mo=8xl0^‘*,  and  the  Tadzhik- 
Sinkiang  1974  event  with  Mo  = 1.7xlo2^. 
On  an  Mq  vs  Mg  diagram,  these  data  tend 
.to  show  higher  apparent  stresses  (nd  > 
100  bars)  for  the  intraplate  events,  as 
Kanamori  and  Anderson  (1975)  have  found 
in  their  data.  The  Brune's  AR  method 
has  been  calibrated  "using  12  events  with 
independent  M estimates  that  have  been 
recorded  on  VWSSN  long  period  stations 
over  predominantly  oceanic  paths.  A 
total  of  15  stations  in  Pacific,  Atlan- 
tic, and  Indian  Oceans,  at  epicentral 
distances  up  to  17000  km  are  used  - at 
least  two  stations  per  event.  An  epi- 
central distance  correction  based  on  a 
more  detailed  Q structure  than  Brune's, 
as  .well  as  an  added  correction  for  in- 
strument response  must  be  applied  to 
give  a nearly  linear  Mq  vs  AR  relation. 
Aftershocks  are  the  largest  source  of 
error  and  give  moment  estimates  with  a 
factor  of  2 error  for  events  with  mo- 
ment S 10^'  dyne-cm  and  a factor  of  3 
error  with  moment  determinations  > 10^^ 
dyne-cm. 
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RAYLEIGH  OCEANIC  GROUP  VELOCITY  DIS- 
PERSION: A PREDICTIVE  MODEL  DERIVED  FROM 

TECTONICS 

Frederick  J.  Mauk  (Department  of  Geology 
and  Mineralogy,  University  of  Michigan, 
Ann  Arbor,  MI  48109) 

A functional  relationship  betv/een  sea 
floor  age  and  group  velocities  of  funda- 
mental mode  Rayleigh  waves  for  periods 
from  30  to  100  seconds  has  been  empiri- 
cally determined.  This  10%  systematic 
increase  in  group  velocity  v;ith  increase 
in  age  is  used  to  define  a bathymetry 
corrected  age-velocity  matrix.  Multipli- 
cation of  this  age-velocity  matrix  and  a 
5*x5®  matrix  encoded  with  sea  floor  age 
has  yielded  Rayleigh  group  velocity  con- 
tour maps  at  periods  of  30,  40,  50,  60, 
70,  80,  90,  and  100  seconds  for  the 
world’s  oceans.  Comparison  betv;een 
predicted  and  observed  group  velocity 
dispersions  using  these  maps  yields 
velocity  residuals  consistently  less 
than  0.02  km/sec. 
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PREDICTION  OF  OCEANIC  RAYJ.EIGH  WAVE 
CROUP  VELOCITIES:  ONE  HIT  AND  SEVERAL 
NEAR  MISSES 

Frederick  J^.  Matik  { University  of  Michigan, 
Department  of  Geology  and  Mineralogy,  Ann 
Arbor,  Michigan  48109) 

A fundamental  mode  oceanic  Rayleigh  wave 
group  velocity  predictor  based  on  sea  floor 
*ge  has  been  described  by  Mavde  (EOS,vol  50, 
no.  10,  1976).  Application  of  this  empirical 
model  to  random  paths  across  the  Pacific 
basins  yield  group  velocity  residuals  for  per- 
iods from  30  to  100  seconds  which  are  con- 
sistently less  than  0.02  km/sec.  Similar  ran- 
dom path  tests  of  the  Pacific-derived  model 
have  beer,  conducted  for  the  North  Atlantic, 
South  Atlantic  and  Indian  Oceans.  Observed 
group  velocities  for  these  three  ocean  basins 
were  consistently  faster,  frequently  greater 
than  0.20  km/sec,  than  predicted  by  the  model. 
The  largest  residuals  occurred  in  the  North 
Atlantic  Ocean,  the  smallest  in  the  Indian 
Ocean.  This  suggests  :hat  models  01  Rayleigh 
group  velocities  based  solely  on  sea  iloor  age 
•re  incomplete.  The  additional  requisitt  model 
correctio.i  to  predict  accurately  the  group  vel- 
ocities fer  all  ocean  basins  appears  to  be  em- 
bodied in  the  spreading  rates.  Since  the  func- 
tional relationship  between  group  velocities 
and  sea  floor  age  appears  to  be  related  to 
the  thernial  evolution  of  the  oceanic  litho- 
sphere, the  observed  \elocity  variations  among 
these  ocean  basins  suggests  that  the  mantle 
heat  flow  into  the  base  of  the  lithosphere  is 
non-uniform. 
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O)  Pa  and  Sa  Phases  Revisited  by  Frederick  J.  Mauk  sub- 
«itted  for  the  47th  Annual  Meeting  of  the  Eastern  Section  of  the 
Seismological  Society  of  America 

"Examination  of  four  months  of  seismograms  from  the 
White  Pine  HGLP  observatory  has  revealed  the  occurrence 
of  vertical  mode  Pa  and  Sa  phases  for  an  azimuthally 
restricted  population  of  teleseismic  events.  These  earth- 
quakes range  in  distance  from  50  to  106  degrees  and  in 
depth  from  30  to  60  kilometers.  The  predominant  periods 
are  15  and  25-50  seconds  for  Pa  and  Sa , respectively.  The 
apparent  velocities  of  both  Pa  and  Sa  show  strong  dependence 
on  the  type  of  lithosphere  through  which  the  waves  pass.  Pa 
ranges  from  8.05  km/sec  to  8.28  km/sec  and  Sa  from  4.58  km/sec 
to  4.54  km/sec.  The  predominant  north-south  azimuths  to  the 
events  yielding  Pa  and  Sa  display  a remarkable  correlation 
with  regions  of  thick  lithosphere  proposed  by  Chapman  and 
Pollack  (1975).  This  suggests  that  the  "whisper  gallery" 
effect  of  Press  and  Ewing  (1955)  rather  than  the  low  velocity 
zone  wave  guide  of  Caloi  (1953)  is  the  propagation  mechanism." 

(2)  A Finite  Element  Model  of  Recent  Turkey  Plate  Movements 
by  Barbara  Radovich  Killiams  presented  at  the  AGU  meeting  in 
Washington,  D.C.,  June,  1975  and  submitted  for  publication  to  the 
Journal  of  Geophysical  Research. 

"Plate  tectonics  in  the  eastern  Mediterranean  shows  the 
African  and  .Arabian  plates  moving  northward  with  respect  to 
Eurasia  and  the  Turkish  and  .Aegean  microplates  moving  rapidly 
in  a direction  perpendicular  to  these  major  plate  motions  in 
a westerly  direction.  Construction  of  a finite  element  model 
of  the  plates  in  this  region,  with  the  microplatcs  loaded  oi\ly 
on  their  edges,  verifies  McKenzie's  hypothesis  that  the  .Arabian 
and  African  plates  colliding  with  Eurasia  are  wedging  the 
Turkish  microplate  westward.  The  motions  of  the  microplates 
in  the  model  are  very  sensitive  to  the  boundary  conditions, 
and  the  Turkisli  part  of  the  model  moves  westward  only  when 
the  Turkey  and  Arabian  plates  arc  uncoupled  from  Eurasia. 

An  easterly  movement  of  the  microplatcs  occurs  when  the 
boundaries  of  the  model  arc  coujilcd  to  Eurasia.  The  smaller 
microplatc  areas  of  the  model  seem  to  act  move  like  the  ideal 
of  a rigid  plate  tiian  the  larger  major  plate  areas.  The 
Border  zone  activity,  on  the  south  border  of  the  Turkish 
microplatc,  is  predicted  to  be  inversely  proportional  to 
the  activity  on  the  Anatolian  Fault.  In  general  the  fit 
to  the  observed  plate  motions  for  the  Turkish  microplate 
arc  good  but  the  model  shows  that  the  Aegean  part  of  tl\c 
model  may  fit  the  observed  seismic  zones  without  really 
moving  to  the  southwest,  but  rotating  counterclockwise." 
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and 

Steven  G.  Henry 


Preface 


The  following  pages  contain  a listing  of  seismic  events 
recorded  at  the  White  Pine  High  Gain  Long  Period  Observatory 
for  the  period  July  1974  through  April  1975.  During  this  period 
the  observatory  was  in  a state  of  instrumental  flux  both  in 
amplifier  type  and  recording  mode.  To  appear  in  the  listing, 
all  events  had  to  exceed  a recorded  signal  to  noise  ratio  of 
2:1  in  Rayleigh  wave  signature.  Only  vertical  seismograms  are 
available  for  this  time  period.  The  exact  arrival  times  of 
body  phases  for  these  events  have  not  been  included  to  encourage 
other  interpreters  to  examine  the  seismograms  more  thoroughly. 
All  of  these  early  grams  may  be  obtained  for  study  by  writing  to 
the  Seismological  Observatory,  Department  of  Geology  and 
Mineralogy,  University  of  Michigan,  Ann  Arbor,  Michigan  48104 
or  calling  (313)  763-8069. 

A brief  description  of  the  station  is  included  in  this 
catalog,  but  data  users  are  encouraged  to  read  the  semi-annual 
Technical  Reports  "Installation  and  Operation  of  a High  Gain 
Long  Period  Seismograph  Station"  ARPA  Order  No.  1827,  Contract 
No.  F44620-73-C-0050  for  details  regarding  alteration  of  the 
system. 


station  Description 


Seismic  Station  WPM  is  comprised  of  a throe  component  high 
gain  (100  K)  long  period  (To  " 30  sec)  seismic  system.  The 
station  is  divided  into  a subsurface  instrument  vault  and  a sur- 
face recording  facility  which  are  linked  via  a hard  wire  teleme- 
try system. 

The  seismometers  have  been  installed  89  meters  below  the 
surface  in  a vault  provided  in  section  31-A  of  the  White  Pine 
Copper  Company  mine,  White  Pine,  Michigan  (see  locator  map  on 
page  4 ).  They  are  sitting  in  specially  constructed  pressure- 
temperature  sealed  tanks  on  the  uppermost  horizon  of  the  Copper 
Harbor  Formation  at  an  elevation  of  193.5  meters  (see  geologic 
column  on  page  5).  All  signals  from  the  seismometers  are  electron- 
ically amplified  and  FM  multiplexed  to  assigned  carrier  frequencies 
to  transmit  them  to  the  surface  recording  site.  Remote  boom  re- 
positioning and  calibration  of  the  instruments  is  accomplished  from 
a subsurface  control  panel  located  in  vault  1 (see  subsurface 
station  schematic  and  floor  plan  on  page  6).  All  recording  of  the 
seismic  signals  is  done  at  the  surface  installation. 

Tests  of  the  subsurface  station  site  environmental  stability 
have  shown;  (1)  it  is  located  in  a portion  of  the  mine  exhibit- 
ing low  cavity  convergence  (40.002  inches/year);  (2)  it  is 
located  in  an  area  having  little  microseismic  activity  but  exhib- 
iting substantial  4 hertz  background  noise  during  active  mining 
periods;  (3)  it  is  not  microbarometr ically  stable  and  exhibits 
a seasonal  fluctuation  in  humidity  ranging  from  99%  in  the  summer 
to  65%  in  the  winter;  (4)  it  is  thermally  stable  to  ^ 1°F  at 
55  F but  does  exhibit  a seasonal  decrease  of  5°F  during  winter 
months. 

The  surface  recording  facility  is  located  in  the  north  end  of 
the  Bechtel  Building  of  the  White  Pine  Copper  Company.  Signals 
received  from  the  subsurface  are  demultiplexed  and  separated  and 


then  recorded  in  both  analog  and  digital  form.  Analog  signals 
are  recorded  directly  on  helicorders.  Simultaneously,  velocity 
and  displacement  analog  information  is  transformed  to  a digital 
format  and  recorded  on  magnetic  tape  with  the  Astrodata  data 
logger  system.  All  digital  tapes  must  be  processed  by  the 
Albuquerque  Seismological  Center  before  they  are  computer  com- 
patible. 

Helicorder  records  are  changed  daily;  magnetic  tape,  bi- 
monthly. All  processing  of  the  data  is  done  at  the  University 
of  Michigan  Seismological  Observatory. 


station  Facilities 


Location 

White  Pine  Copper  Company 
White  Pine,  Michigan  49971 

Station  Abbreviation 
WPM 

Station  Director 

Dr.  James  T.  Wilson 

Station  Address 

c/o  Fred  Mauk 

Seismological  Observatory 

Department  of  Geology  and  Mineralogy 

C.C.  Little  Building 

University  of  Michigan 

Ann  Arbor,  Michigan  48104 

Telephone 

Mine  (906)885-5111,  ext.  321 
UofM  (313)763-4069 

Geographic  Coordinates 

45°45'05.5"  North  latitude 
89°33'17.75'’  West  longitude 

Station  Elevation 
193.5  meters 

Depth  of  Sensors 
89  meters 

History  of  Station 

The  White  Pine  high  gain  long  period  seismograph  station 
was  established  in  the  White  Pine  Copper  Company  mine  during 
1973  and  1974  with  funds  provided  by  the  Advanced  Research 
Projects  Agency  under  Contract  No,  F44620-73-C-0060.  First 
installed  as  an  experimental  station  in  1973,  WPM  became  the 
eleventh  HGLP  observatory  in  1974.  For  fvirther  installation 
information,  refer  to  ARPA  Technical  Reports  Contract  No. 
F44620-73-C-0060. 
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All  events  for  the  period  July  1974  through  November  1975 
were  recorded  with  a Geotech  7505A  long  period  seismometer 
(To=30  seconds) . The  signal  was  amplified  using  a Geotech 
AS-650  solid  state  long  period  amplifier  and  multiplexed  for 
telemetry  purposes.  The  signal  was  demultiplexed  using  a 
Geotech  XD-410  discriminator  and  recorded  at  two  gain  levels 
on  Esterline  Angus  S-601-S  Strip  Chart  recorders.  The  gains  are 
split  fay  one  order  of  magnitude  8k  and  88k.  The  timing  system 
is  a Springnether  TS-100.  All  events  from  December  1974  through 
April  1975  were  recorded  on  a Geotech  RV-301-B  rectilinear 
helicorder  at  a gain  of  5k.  '~he  timing  system  is  a Geotech 
TG-120.  For  further  information  refer  to  Technical  Reports 
ARPA  Order  NO,  1827  contract  No.  F4462D-73-C-0060 . 
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Az  and  Baz  measured  CW  from  North 
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Az  and  Baz  measured  CW  from  North 


The  following  pages  are  plots  of  events  recorded  at  the 
White  Pine  Observatory  on  a monthly  basis.  Each  plot  Is  an 
azimuthal  equidistant  projection  centered  on  White  Pine. 

Each  tic  mark  represents  a 10®  distance  increment  with  the 
outer  circle  at  150®.  The  size  of  the  "X”  plotted  is  de- 
pendent on  the  reported  mb  magnitude  of  the  National  Earth- 
quake Information  Service.  When  several  events  were  record- 
ed from  the  same  locality,  all  subsequent  events  were  dis- 
placed slightly  to  make  them  visible. 


-26- 


NPM 


MRRCH  EVENTS 


1 


UNCLASSIFIED 


^ TITLE  (and  Subtitle) 

n j ^IGH  £AIN  j^NG  PERIOD  SE 


OP  This  page  flVhen  Date  Entered) 


DOCUMENTATION  PAGE 


(2.  GOVT  ACCESSION  NO 


SEISlVDGRAPH  STATION^ 


ij  Final  Apr-  d 

I— 


lywilson  ^ 


James  "(^/Wilson 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Institute  of  Science  and  Technology 
University  of  Michigan 
Ann  Arbor.  MI  481 C5 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

hRPA/iJMR 

140C  Uilson  Joulevard 
Arlington,  Va  22209 


READ  INSTRUCTIO.NS 
fJEFOKE  COMPLETI.N’G  FORM 


3.  RECIPIEUT'3  CATALOG  NUMBER 


O*" 


cor'^"i"irT  n°  ' ~ I J 

FA462)/)-73-C-f6l^6g(^f^'^/j h 

CU 


10.  PROGRAM  element, 

AREA  & WORK  UNIT 

AO  3291-1  62701 E 

7F10 


14.  MONITORING  AGENCY  NAME  d ADORESS^^// d/7/erenf  Irom  Corrtrolling  Office) 

AF0SR/i<P 

Bolling  AFB,  Bldg. #410 
Wasn.,  O.C.  20332 


13.  NUMBER  OF  pZs«  ' -* 


104 


15.  SECURITY  CLASS,  (of  this  report) 

Unclassified 


15«.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  fof  Ihis  Report; 


Approved  for  public  release,  d.'.stribution  unlimited 


17.  OISTRiauTlON  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 


18.  supplementary  NOTES 


19.  KEY  WORDS  fCont^nue  on  reverse  side  if  necessary  and  <denf//y  by  block  number) 


ABSTRACT  (Continue  on  reverse  side  if  necessary  and  identify  bv  block  number) 

This  final  report  on  Contract  No.  F44620-73-C'-OC60,  April  1,  1973  - 
March  15,  1978,  describes  the  work  carried  out  under  the  terms  of  the  contract 
A High  Gain  Long  Period  (HGLP)  seismograph  station  was  installed  in  the  White 
Pine  Copper  Mine  in  Ontonagan  County,  Michigan.  The  station  went  through  a 
series  of  stages  starting  with  a long  period  vertical  recording  through  a 
phototube  amplifier  on  a strip  chart  and  ending  as  an  essentially  standard 
HGLP  station.  Because  of  the  difficulties  of  servicing  and  recording  withDut/! 

J.gi  Ufi ill 


DO i473r  EDITION  OF  1 NOV  65  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  fUTien  Dttm  Enf*r«d; 


J 


ili: 


SECURITY  CLASSIFICATION  OF  THIS  PACECHTian  Data  Bnltrmd) 


telemetry,  the  station  was  closed  at  the  end  of  1976.  Much  useful 
information  regarding  the  engineering  and  installation  of  HGLP  stations 
in  mines  was  obtained.  The  station  produced  many  months  of  useful  long 
period  data.  Concurrently  with  the  station  installation  two  extensive 
research  efforts  in  long  period  seismology  were  carried  on.  One 
developed  a Rayleigh  wave  velocity  model  for  ocean  basins  and  the  other 
applied  long  period  P wave  spectra  to  the  determination  of  seismic  moments. 


'■  > . .A-  r 


SECURITY  Cc  ASSiriCATlOR  OE  Tw*  f Ertmrmd) 


